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FOREWORD 


This final report documents the technical studies conducted by Ford Aerospace 
& Communications Corporation, Aeronutronic Division under Contract 955115 to 
the California Institute of Technology Jet Propulsion Laboratory (JPL) in 
Pasadena, California. The JPL Technical Manager was Mr. J. R. Womack. 

This is a three volume report prepared by the Aeronutornic Division. Subcon- 
tractors were the WDL Division of Ford Aerospace & Communications Corporation, 
Palo Alto, California; United Stirling of Sweden (USS) , Malmo, Sweden; Sund- 
strand Energy Systems, Rockford, Illinois. 

The WDL Division was responsible for the concentrator and electrical subsystems. 
USS provided information on Stirling engines, and Sundstrand supplied infor- 
mation on organic Rankine-Cycle Engines. Additional supporting information 
was provided by Garrett AiResearch Manufacturing Company, Phoenix, Arizona 
(closed-cycle Brayton engines); Solar Turbines International, San Diego, 
California (open-cycle Brayton engines); and Williams Research, Walled Lake, 
Michigan (open-cycle Brayton engines). Also, the following divisions of the 
Ford Motor Company provided expertise: Glass Division, Scientific Research 

Laboratory, and the Manufacturing Planning Group. 


The key personnel for the studies documented in this final report are listed 
below: 

• Aeronutronic Division, Ford Aerospace & Communications 
Corporation 

N. L. Cowden - Program Manager 

R. L. Pons - Technical Manager and Systems Analysis 
T. B. Clark - Optics and Assistant Technical Manager 

D. B. Osborn - Thermodynamics 

E. D. Avetta - Design 

D. C. Jackson - Structures 

• WDL Division, Ford Aerospace & Communications Corporation 

H. J. Sund - Subcontract Manager 

I. E. Lewis - Concentrator Subsystem 

J. L. Knorpp - Electrical Subsystem 

• United Stirling (USS) 

W. Percival - Consultant to Aeronutronic and Marketing 
Representative of USS 

Y. Haland - Program Manager for Solar/Stirling 
Applications 








Sundstrand 


M. Santucci - Principal Investigator , Organic Rankine 
Engine 

• Garrett 

L. Six - Principal Investigator, Closed-Cycle Bray ton 

Engine 

• Solar Turbines International 

M. Gramlich 

• Williams Research 

R. Mandel 

• Ford Motor ComDanv Manufacturing Planning Group 


W. Nagle 


T. B. Clark was the editor of these reports. 
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SUMMARY 


This report summarizes the results of meetings held with control equipment 
suppliers (3asler Corp.), the local electric utility (Southern California 
Edison) and the engine manufacturer (United Stirling of Sweden) , to establish 
the basic control configuration for the baseline solar SPS. An overall control 
arrangement is established and various interface requirements are briefly 
considered. 
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BACKGROUND 


The objectives of Che short study reported on in this TR were to examine 
1) the special requirements for control of a Stirling cycle engine when 
employed in a solar application, and 2) the requirements for interfacing a 
solar-powered eugine with an electrical utility grid. 

A fuel-powered Stirling engine, as developed by United Stirling- Sweden (USS) 
for automotive applications, incorporates two separate loops for control of 
output power (torque and speed) and efficiency. Figure 1 shows the power 
control loop; engine shaft output is controlled by controlling the pressure 
and hence the mass of the working fluid* within the engine. It is known as 
mean pressure level control. The system consists of a high pressure reser- 
voir, a compressor, a servo-driven valve** and associated plumbing that per- 
mits the working fluid to be moved into or out of the engine at will. Note 
that this technique of "inventory" control is characteristic of closed cycle 
heat engines. 

Figure 2 shows the second control loop, which provides for metering both air 
and fuel flow (modified Bosch K-Jetronic system) to the burner head of the 
engine in order to keep head temperature (actually tube wall temperature) 
constant -- and hence maintain constant engine efficiency — independent of 
the operation of the power control loop. 

SPECIAL SOLAR REQUIREMENTS 

It is both impractical and undesirable to try to vary the sun's input to the 
engine -- as the automotive Stirling varies the amount of fuel burned. The 
impracticality stems from the difficulty of either shuttering the solar input 
or of bypassing part of the sodium vapor, e.g., to a condenser and back to the 
receiver/boiler unit. The undesirability stems from the reduced power output 
associated with this method of control. It seems obvious, therefore, that the 
control scheme must be configured to accept all of the solar energy directed 
to the receiver. One way to accomplish this is to "slave" the existing power 
control loop to engine head (or receiver) temperature; this arrangement permits 
output power to "follow" the solar input power curve while maintaining constant 
system temperature. There are special requirements for system start-up and 
subsequent grid synchronization, however, as discussed in the following para- 
graphs. The information presented in these paragraphs is based on extensive 
discussions with engineering personnel of the Basler Corp. (a major manufac- 
turer and installer of power plant control equipment) and the Southern Cali- 
fornia Edison Company. Preliminary discussions were also held with engineering 
personnel representing the engine manufacturer, USS. 


* For a fuel-powered engine, the working fluid is Hydrogen; for the solar 
application, it is Helium. 

** The power control valve is a hydraulic slide valve unit designed and built 
by USS; it is driven by a Moog Series 76 electro-hydraulic servo valve. 
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OVERALL POWER SYSTEM 


The block diagram for the overall system is shown typically in Figure 3. 

Each solar module (collector, Stirling engine and alternator assembly) will 
have separate automatic and manual override controls to place each alternator 
''on-line" by actuating a unit power contactor. For maintenance, repair, etc., 
each set of modules could also have a power contactor (both automatic and 
manual modes) to place the entire string, e.g., line "A" of Figure 3, on the 
grid. (The line power contactors are an optional feature.) Also, the energy 
storage system will be controlled to supply stored energy per overall system 
requirements. The microprocessor will provide automatic control for all 
modules except for the automatic "on-line" and "off-line" control equipment 
built into each module. Energy delivery sensing and recording equipment will 
be supplied by the utility company. The master control and transfer unit 
will include a mimic bus and status indicators. All alternators will operate 
in the "scramble" mode (each alternator independently "scrambles" to g-'f onto 
the utility grid as soon as it can) and each alternator voltage regulator will 
operate with parallel droop compensation. All alternators are delivering 
power to an "infinite" bus, i.e., the voltage and frequency are independent of 
the solar system power level. 

STIRLING ENGINE OPERATION 

As mentioned above, it is necessary that the engine remove power from the 
receiver at essentially the same rate as the net solar input to the receiver. 
This automatically assures constant temperature operation, which in turn 
provides high overall system efficiency and minimizes the deleterious thermo- 
structural effects of thermal cycling. The only direct means for controlling 
engine output is through the existing power control subsystem. This is 
accomplished by driving the power control valve so as to satisfy the require- 
ments for different modes of operation, e.g., 1) start-up, 2) grid synchroni- 
zation, 3) on-line operation and 4) emergency shutdown. A preliminary approach 
to establishing the control sequences and associated hardware is given below. 
The approach is predicated on the following observations: 

• The combination of receiver and heat transport components comprise 
a very slow response system (typically 1 to 2 minutes to show 
significant temperature change due to large step changes in either 
input solar flux or in power removed by the engine) 

• The engine/alternator combination is a comparatively fast response 
system (idle to full power in less than 1 second) 

• It is undesirable to try to pull high power from the engine at 
head temperatures less than about 100°C-200°C below steady-state 
operating temperature. In a closed system as employed the rapid 
drop in sodium vapor pressure and density with lowered temperature 
results in flow choking in the vapor pipe connecting the receiver/ 
boiler to the head of the engine. An unreasonably large increase 
in pipe diameter would be required to avoid choking at the vapor 
mass flow rates rresponding to operation near rated power. 
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« For a system with a large number of modules, it is desirable to get 
all engines up to speed and quickly locked onto the utility grid 
at low power level, then bring each engine up to full power as soon 
as the operating temperature range is reached. 

Careful consideration of these points suggests the following sequence: 

1. Start each engine as soon as there is sufficient temperature and 
input solar flux for self-sustained operation. The engine will 
operate (unloaded) at idle speed, developing sufficient power to 
overcome internal work (friction and compression) as well as the 
drag of the alternator. 

2. Move the power control valve to increase Helium pressure in the 
engine, thereby increasing engine torque and speed until alternator 
output frequency, voltage and phase are matched to grid values. 
Close line switch (contactor) and lock onto grid at lower power 
level. The grid will hold the engine at essentially constant 
speed. 

3. When engine head temperature reaches the proper value, move the 
power control valve further and increase engine torque (speed is 
constant) until the energy withdrawn from the receiver matches the 
net input from the sun. As the solar power input increases or 
decreases, modulate the power control valve to maintain constant 
temperature (within the appropriate control band) . 

4. In che event of an emergency, e.g., sudden unloading of the engine, 
overspeed damage would be avoided by quickly moving the power 
control valve so as to release the working fluid from the engine 

to the reservoir, thereby shutting the engine down. 

Implementation of this sequence of operation is proposed as follows: 

Start-up 

At shut-down on the previous day, the remotely- actuated blocking 
valve was closed to prevent heat leak to the engine overnight via the 
sodium vapor line. Also, most of the Helium was evacuated from the 
engine to minimize seal leakage overnight. When the sun reaches the 
minimum elevation angle achievable with the collector, identifiable 
either by a sun-tracking signal or by signal from a programmable 
timer, 1) the module is moved into position, 2) the cover of the 
cavity receiver is rotated into the stowed position (letting the 
focussed solar beam into the aperture) , 3) the power control valve is 
moved to permit engine pressurization with Helium to 30-50 atm, and 
4) the blocking valve is opened to permit the still-hot* sodium vapor 
to warm up the vapor pipe and the engine heater head assembly. 


*Normal nighttime conditions serve to keep the sodium temperature in the 
receiver well above 250°C; for extended bad weather, however, the sodium 
can freeze and the engine should not be started until the sodium is melted 
and raised to at least 250 3 C. 
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Upon receipt of sensor signals that 1) verify that head temperature is 
> 250’C, 2) verify that the sun is putting at least "threshold" energy into 
the receiver and 3) verify engine Helium pressure level, the starter motor 
is engaged. 

With the power control valve set in the idle power position, the engine will 
come up to 750 rpm idle speed in approximately 6 seconds and the starter 
motor will automatically disengage. Engine speed will hold constant at 
750 rpm by action of a shaft governor which modulates the power control valve; 
it is part of the existing power control system. The engine is now ready to 
be placed on-line while the sun continues to warm up the receiver. 

Synchronization 

When the head temperature reaches a pre-selected value (~ 500°C), the power 
control valve is placed under command of the grid synchronizer (typically a 
Basler PRS-370 and associated equipment) which moves the valve to increase 
Helium pressure within the engine, causing a simultaneous increase in engine 
shaft torque and speed (increased power). Once the correct speed is reached, 
the synchronizer "dithers" the valve (and thus engine speed) to obtain cor- 
rect phase; upon verification of voltage, frequency and phase match with the 
grid, the synchronizer activates a unit power contactor and the system is 
locked onto the grid, with the engine speed held essentially constant and 
the power valve in a "minimum" power setting. With this action, the synchro- 
nizer releases control of the power valve to the engine head temperature 
sensor, as shown in Figure 4, and the system is now ready for normal on-line 
operation*. 

On-Line Operation 

Engine heater head temperature will continue to rise since the net solar 
input power is much greater than the power removed by the engine. The engine 
power setting will remain at the "minimum" position until the temperature 
rises to a pre-determined value, e.g., above the steady state operating point 
but within the temperature deadband. At this point the temperature controller 
moves the power valve to bring the engine up to a torque level such that more 
energy is removed from the receiver than is supplied by the sun and the tem- 
perature will begin to drop. From this point on the power valve is modulated 
to hold temperatures within the control band as the solar flux varies during 
the day. Figure 5 shows temperature, speed and torque histories for a typical 
starting sequence. 

Emergency Shutdown 

As indicated previously, an overspeed sensor is keyed to the engine drive 
shaft and switches the engine power control valve into a "dump" position so 
that the working fluid is rapidly removed from the engine and the power is 


*Further FACC power control subsystem analyses suggest that an alternate 
approach employing the central microprocessor for sequential synchronization 
is less costly. The actual synchronization procedure, however, is unchanged. 
Documentation of the results of these analyses is forthcoming. 
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reduced to idle level (or zero) to prevent damage in the event of 
inadvertent, sudden loss of load from the grid. This is similar to 
the "short-circuit" feature currently incorporated in the automotive 
power control. Other events would occur simultaneously with the 
power dump operation, e.g., the concentrator azimuth/ elevation drive 
system would be signalled to defocus the complete module and the unit 
power contactor would open to make certain that the alternator is 
off-line. 

United Stirling has indicated that it may be desirable to modify the existing 
"short-circuit" technique if engine shutdown in ~ 20 seconds is acceptable. 
Furthermore, since the solar operation does not have the numerous start/stop 
characteristics of the automotive operation, some power can be saved by adding 
a small low-pressure reservoir on the suction side of the compressor, and 
running the compressor only intermittently, i.e., just to keep the high 
pressure reservoir charged up (typical of air brake operation on trucks). 

Engine load is supplied by the alternator, whose operation is described in 
the following paragraphs. 

ALTERNATOR 

The operation of the alternator is considered in terms of its operating modes, 
control equipment/ operation and its interfaces, as described below: 

1. Basic Operation 

The alternator is directly connected* to th?-. Stirling engine; a shaft 
coupling will be provided for assembly/disassembly as required. The power 
input to the alternator is thus the power output of the Stirling engine 
which is proportional to the product of Torque X Speed. As previously stated, 
once the alternator is connected to the utility grid its speed is held fixed 
by the grid since speed is proportional to frequency. Thus, as the engine 
power output is varied by the control system to accommodate varying solar flux, 
only the shaft torque will vary. The alternator power output is in turn pro- 
portional to the product of Voltage X Current since its output voltage is 
held fixed by the utility grid voltage, the current will vary in proportion to 
the input torque variation. 

2 . Controls 

Alternator/ engine control arrangements are shown in Figure 6. A number 
of synchronizer circuits are available to synchronize the alternator with the 
electric utility grid before connecting to it. The Basler PRS 370 will pro- 
vide the proper phase matching functions between the alternator and the line 
by controlling the speed of the Stirling engine, as discussed previously; its 


*System studies indicate optimum engine speed is 1800 rpm, which is the 
selected value for the alternator; thus, no gearbox is necessary. 
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approximate cost is $520. The Basler PRS 150 Sync Check, at an additional 
cost of $918, mav be desirable. It will not allow the alternator on-line 
unless both phase and voltage are within close tolerances. 

Southern California Edison requires Volt-Amperes Reactive (VAR) control to 
minimize VAR losses in the electrical system. The alternator and the load 
usually have some inductive reactances which result in load conditions 
having both in-phase power and reactive power as shown in Figure 7 below: 



The horizontal leg of the triangle represents true electrical power (KW) 
delivered; the vertical leg is proportional to the kilovolt-ampere product 
and represents reactive power. The cosine of the angle, 9, is the power 
factor; to prevent excessive VAR losses, the power factor must be kept above 
approximately 0.8. When the alternator is connected to the utility grid, 
this is accomplished by varying the alternator field to control the KVAR/KW 
ratio. Basler makes VAR control circuits for most alternators. They can 
supply an existing VAR control as well as a voltage regulator to go in the 
alternator circuit as shown previously in Figure 6. 

A reverse power control circuit is also required to prevent utility grid 
power from entering the alternator. A Basler BE-1-32 unit (or equivalent) 
must be included. 

3 . Multiple Alternator Operation 

As discussed previously, it is desirable to connect the alternator to 
the grid as quickly as possible to provide load for the engine. By operating 
all alternators in the scramble mode, each alternator will go "on-line" as 
soon as it is rotating at the required speed (frequency) and phase. This is 
the normal mode for utility operation where the alternators are being operated 
into an "infinite" bus. For example, at the San Onofre Nuclear Plant, 12 
steam turbine-electric generating units -- at 50 MW each -- are "scrambled" 
onto the SCE/SDGE grid, without introducing dynamic instability in the grid. 

The voltage regulators will be connected in the parallel droop compensation 
mode. In this mode, the regulator circuits have interconnections between 
all alternators. Without these circuits, one or more alternators might have 
excessive currents in the windings resulting in an overload condition. 
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4. Ucilitv Interface 


The utility interface will be as indicated in Figure 8. The utility will 
provide sensing circuits as indicated and the power contactors will be opened 
ac any abnormal operation. SCE indicates that there may be no warning before 
they disconnect, so it is essential that we be able to de- focus the solar col- 
lectors sufficiently fast to prevent damaging temperature overshoot in the 
receiver and in the engine heater head. 

The solar power system will need lightning protection; SCE would provide this 
protection in their package, if necessary; further analysis of this issue is 
required. 

In consideration of Phase II of the program, SCE indicated that SPS power -- 
from one or more test units -- could be worked into the Aeronutronic plant 
electrical grid if desired. The special interface equipment would cost 
roughly $200 per month. 

CONT U S IO NS 

A preliminary analysis of the SPS control requirements has been made and an 
initial system mechanization identified. More detailed analysis is required 
to definitize the system. For example, the following information is required: 

1. Transfer characteristics of the current Stirling engine power 
control valve 

2. Engine head temperature sensor characteristics 

3. Data requirements for the microprocessor 

4. Alternator characteristics 

5. Other interface data, e.g., the energy storage system, the master 
control and transfer unit (including the mimic bus and status 
indicators) , the energy indicator and recorder, the energy 
transfer unit 

Also, additional analysis should be carried out to address the control require- 
ments for accommodating engine operation during transient cloud passage. 
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SUMMARY 


The Garrett CCPS-40 Gas Turbine is mathematically modeled to permit deter- 
mination of annualized SPS power output and to identify working fluid mass 
flow rates required for sizing the solar receiver. 
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BACKGROUND 


This report summarizes results of a preliminary analysis carried out for 
a closed Brayton cycle engine -- specifically the Garrett CCPS-40-1 Gas Turbine 
shown in Figure 1 -- as an initial step in evaluating its applicability to the 
solar SPS concept. The basic objectives of the analysis were 1) to determine 
engine output (shaft) power, P q , as a function of input power, P^ and ambient 
temperature, T , and 2) to establish working fluid mass flow requirements. 
Objective (1) was formulated to permit determination of annualized SPS power 
output; objective (2) was formulated to permit sizing of the solar receiver. 

Selection of a closed cycle engine for initial study was predicated on a 
recommendation by Garrett personnel for the following reasons: 

1) The engine can be operated at constant speed and constant turbine 
inlet temperature (TIT) over the entire power range. Power variation 
is achieved by inventory control, i. e. by varying the amount of working 
fluid in the engine. Engine efficiency is nearly constant over the 
power range. Without an expensive inf initely-variable ratio gear 

box, an open-cycle engine would require TIT variation for power con- 
trol at constant speed and part-load efficiency would be substantially 
reduced. 

2) By operating at high working fluid pressure, effectiveness of the 
heat exchangers is improved, resulting in smaller, lighter, less 
costly units than would be required by an open cycle machine. 


APPROACH 

The approach taken here employs well-established theoretical cycle analysis 
to generate an equation relating P^, P Q and T a . The cycle analysis, however, 
is modified to accommodate component performance data derived from tests con- 
ducted by the U. S. Navy on the CCPS-40 engine as reported in Reference 1. 

Performance estimates are also made for the CCPS-40 engine with future 
improved component efficiencies as suggested by Garrett. 


ANALYSIS 


As derived in the Appendix, a generalized efficiency equation for the 
CCPS-40 engine is 
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■ Recuperator effectiveness (high pressure side) 
e 2 ■ Cooler effectiveness (working fluid side) 


Note that Equation (1) Includes an allowance for bleeding some of the working 
fluid (Argon) off of the compressor discharge in order to cool the bearings, as 
reported in Reference 1. Also, the factor is introduced to accommodate the 
influence of mechanical/thermal losses and Argon leakage; comparison with the 
test data of Reference 1 indicates that ijr * 0.93 gives good data correlation. 

• Compressor Pressure Ratio 

Figure 2 shows engine efficiency as a function of compressor pressure 
ratio, R c , for the CCPS-40 engine as it is presently configured, with the 
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following parameters and component efficiencies corresponding to a nominal 
power output of 25 KW: 


c 

- 0.757 

9 - 

0.935 

\ 

- 0.872 

TIT 

- 816°C (1500 F) 

€ 1 

» 0. 895 

T 

a 

- 44. 6°C (112 F)* 

*9 

* 0. 960 




Figure 2 shows that maximum efficiency occurs in the vicinity of R c ■ 2 ; 
according to Reference 1 the nominal R c for the CCPS-40 engine is 1. 89. Subse- 
quent performance calculations have all been carried out at this value, assumed 
invariant with changes in working fluid pressure/mass flow. 


• Effect of Component Performance 

Reference 1 provides test data on compressor and turbine efficiencies as 
well as on recuperator and cooler effectiveness. The smoothed data for these 
parameters are given in Figure 3, corresponding to the following curve fits: 


% 

- 0.77 - .0005 P Q 

(3) 

It 

- 0.0683 + 0.6941 (P 0 )* 046 

(4) 

*1 

- 0.9249 - 7.917 (10)‘ 5 (P Q ) 1,79 

(5) 

*2 

- 0.98015 - 1. 3474 (10)“ 9 (P ) 5 * 06 

o 

(6) 


Substitution of Equations (3) - (6) into Equation (1) yields engine efficiency 
vs. output power as shown in Figure 4. The comparison with the test data of 
Reference 1 is well within + 107, over the entire power range. The data are re- 
plotted in Figure 5 as P Q vs. P^; Equation (1) is seen to correlate the test 
data quite well. 


• Generalized Performance 

It is assumed that TIT will not exceed 816°C (1500°F) for the foreseeable 
future, primarily due to materials limitations; for R c ■ 1. 89, 0 * 0.935 and 
5> c /m t ■ 1.033, Equation (1) for CCPS-40 engine efficiency becomes: 



where 

T t “ 1 - H t 5 t “ 1 - (0. 204)T) t 

T c “ 1 + " 1+0 * 291/ \ 

Equation (7) is plotted in Figure 6 for ambient temperatures of 10°C, 30°C, and 
50°C. A curve fit to these data yields: 


0.93 


( 7 ) 


P - [0.30375 - . 001275(T )1 (P, + 6) -1 

o a ■* l 

where 



P , P. ~ KW 
o* i 

Equation 8 thus provides the desired relationship between P., T and P . 

X & O 


• Hass Flow Rate 

The power generated by the turbine, P , is: 


C P Vt*’ 


K-cal 

sec 


where 

m t Kg/sec 

C ~ 0. 124 K-cal/Kg°K 
P 



and the power consumed by the compressor, P , is 


( 8 ) 


(9) 
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t 


K-cal 

sec 


( 10 ) 


P 

c 


ra C T ' i 
c p l a c 


The net power output is 


thus P 

o 


P or: 
c 


P 

o 


(4.188) m C T 

C p 3 




( 11 ) 


Equation (11) permits solution of the mass flow rate, m^, as a function of P q 
and T , for the proscribed T, ■ 816°C, R a 1.89, 9 * 0.935 and the curve 

3L C 

fits given by Equations (3) - (6). Typical results are shown in Figure 7, mass 
flow rate is nearly linear with output power. 


POTENTIAL PERFORMANCE IMPROVEMENTS 

If TIT is held to 816°C for the foreseeable future, improvements in engine 
performance must come from improvements in the various components. The "schedule" 
for such improvements is obviously speculative since there are significant cost 
implications. We assume, however, that for the near future, i. e. ca 1982, the 
most expeditious thing to do is 1) to improve recuperator effectiveness, ei, by 
increasing heat exchanger surface area and 2) to increase 9 by increasing flow 
passage size wherever possible, particularly in the interconnecting ducts. 

Both types of improvements obviously require engine repackaging. The performance 
gains corresponding to estimated values of ■ 0.93 and 9 » 0. 95 are shown in 
Figure 8 for a nominal power output of 25 KW. Peak efficiency is increased 
from 25-1% to 28.4%. 

For the far term, ca 1990, the rotating group could be redesigned; the 
additional performance gain achieved by going to estimated values of T| c ■ 0.85 
and ■ 0.92 are also shown in Figure 8; max Tig = 36.77®. 


CONCLUSIONS /RECOMMENDATIONS 


A simple expression has been derived which will permit determination of 
the annual power output of a solar-driven closed Brayton cycle engine, specific- 
ally the Garrett CCPS-40 engine. Working fluid mass flow rates have also been 
determined to permit sizing the solar receiver. 

Estimates of improvements in component performance with further development 
have also been made and the effect on engine performance has been evaluated at 
a single power point. It is recommended that performance expressions similar to 
the aforementioned one be generated over the complete power range to ascertain 
the annualized performance change with the estimated component improvements. 
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It is also recommended that similar analyses be carried out for the open 
Brayton cycle engine before final jvdgements are made regarding the optimum 
SPS engine selection. 


REFERENCE 

1) Duvall, Gordon D. , "Operational Evaluation of a Closed Brayton Cycle 
Laboratory Engine", 11th Intersociety Energy Conversion Engineering 
Conference, 1976. 
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SO 
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• SOLAR INPUT, Q go 

• « / t 3 \ 

Qso “ m t c p^ T 4 ’ t 3^ " m t c p T 4 V 1 " x7 . 


t 3 - t 2 
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*1 = 

T 5 - t 2 

SO 
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T 5 (1 ' 
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T 3 /T 4 

t 2 

- t 7 a ’ 4 i 
4 
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c 
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(A2) 


(A3) 


(A4) 


1 


— — - 


± 




Substitution of (A4) into (A3) gives 




• T T 

n t c p x 4 _ 1 " 'c ( tI T 4 • “ 8 l T t J 


T,/T must now be determined in order to evaluate Equations (Al) and (A5) 

i cl 

only in terns of R c , T a and T 4 : 


T 6 > Tx Ti T 6 

S 2 “ T ' "- — » S0 t" = Ti (1 ‘ S 2> + 6 

A 6 1 a A a a 


OR F 
A a 


/ T 6 N f t 4 \ 

l -2> {-.) ( tT ) T t + 


Fran heat balance on Che recuperator, 

Tg T 3 T 2 

I 5 - I 6 - Tj - T 2 . SO - - i-- 5 +~ 5 



Substitution of Equation (A4) into (A7) yields: 



Substitution of Equation (A 8 ) into (A 6 ) yields: 



WHERE 


Equation (A9) thus becomes: 



• CYCLE EFFICIENCY, 

P fc P 

« • t ~ c 

'CY Q 



WHERE 

I. - 1 - OR.)' k 



Equations (All) and (A10) now fully define the cycle efficiency in terms of 

the variables T , T. , c, , e- f 3 and R . Engine efficiency can be written as 
a 4 Jl i c 

^E " ^CY ^ 


where •) is a correction factor (~ 0. 93 for the CCPS-40 engine) to account for 
mechanical /thermal losses and the effect of working fluid leakage. 


PART 2 



l 

OPEN CYCLE PERFORMANCE MODEL 

ENGINE SCHEMATIC 

The engine analyzed herein is shown in the following sketch. 



• HEAT BALANCE FOR THE RECUPERATOR 
m t (T 3 - T 2 ) - (m t - n^) (Tj - T fi ) 

T 3 - T 2 - (1 - 5)(T 5 - T fi ) 


where 


6 - ctM 

*b 

a * — = By-pass ratio 


m 

M ■ — = Bleed ratio 
in 
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• ENGINE EFFICIENCY, • £ 



E (^) U " Vc (1 " 5)1 ' T c [1 ‘ e c (1 ' 5)1 


(B5) 


(B6) 


where 

e * e = fn(a) (See Recuperator Effectiveness Derivation below) 
c H 


Thus 



(B7) 


Substituting appropriate constants into equation (B7) give the following 
expressions : 


/P V 

/ p i\ 


/ P i\ 

(if) 

- 0.00412T! 

for 

(0 


(B8) 


where 


P = Specified output power, kW 
°r 


P. » — - where “n - ■ 0.28123 

^ C E r 

r 


■ Ambient temperature, K 
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and 




(0.00412)^ 


where 



=■ 1.1781 - 0.0005604(1^) 


* 2.7214 - 0.0054145T 1 



(B9) 


Equations (B8) and (B9) have been used to model the open-cycle Brayton engine 
using program SPEE1 with the Barstow solar data. 

• RECUPERATOR EFFECTIVENESS DERIVATION 
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APPENDIX G 


DETERMINATION OF STATION- KEEPING POWER 


All system calculations carried out to date include an allowance for the 
parasitic power required to operate the SPS. The station power requirement 
is comprised of 1) the power needed to operate each collector and 2) the 
power to operate and environmentally condition the control building. These 
loads were evaluated for both the no- storage case and the case of battery 
storage; the analysis is incorporated in TR-E-4 which is part of this 
Appendix. The results of TR-E-4 have been modified to accommodate more 
recent system design changes. For example, the original storage energies 
of 648 and 7848 kWh no longer correspond to the energies associated with 
annualized capacity factors (ACF) of 0.4 and 0.7, respectively. The data of 
TR-E-4 have been generalized to scale to any level of storage energy re- 
quired, as described below. 

1. NO STORAGE 

The average daily collector load is given in TR-E-4 as 14.85 kWh but this 
value includes power for engine control as well as the engine starting 
battery charger. If these values are removed (they are already incorporated 
in the quoted USS P-75 engine performance) and an average operating day of 
~8 hours employed instead of the quoted 10.14 hours, total collector power 
(at the grid) can be represented by: 

P C0LL ** 1,5 N °* kw 

where N 0 is the number of collectors. 

The building load — which is primarily for air conditioning the 680 ft2 
structure — is given as 125.27 kWh /day in TR-E-4, but this conservative 
analysis was carried out for summer heating conditions. For yearly average 
conditions and an approximate 8 hour operating day, the building power re- 
quirement can be expressed more accurately as 

P BLDG " 8 kW < 2 > 

The station power requirement without storage is thus: 

P 0 (STA) - 1.5 N fl + 8, kW (3) 

Equation (3) is incorporated in the computation of system output based on 
the 15 minute Barstow site tapes. 

2. BATTERY STORAGE 

The station power must be increased to accommodate the increased building 
loads associated with the addition of batteries. The added station power is: 


G-l 


(4) 







E 1 " Wyr 

'./here 

* Added daily energy co building. 

* Constant derived from data of TR-E-4, = 0.177 (1/hr). 

E, = Average daily energy delivered to the grid from batteries, 
kWh. 


- Number of operating days per year. 

Figure G-l presents the additional building energy required for battery 
storage as a function of the battery energy delivered to the grid. The 
summer operation data from TR-E-4 was used to size the air conditioning 
units. This ensures that proper environmental conditions exist within the 
storage building during hot days. However, these data are not appropriate 
for cooler days since the air conditioning equipment will be operated much 
less frequently (or not at all). Thus, the data of TR-E-4 have been ad- 
justed to reflect the average annual energy requirements for the building 
(dashed line in Figure G-l). Approximately 17.7 percent of the energy 
delivered from battery storage to the grid is required to condition the 
storage building on an average annual basis. This percentage corresponds to 
the constant C given in Equation (4). 


► 



r 

\ 
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ADDITIONAL BUILDING ENERGY REQUIRED, E kWh 
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TO 


SUMMARY 


Station auxiliary device loads have been analyzed and found to fall in 
two classes; i.e., those directly associated with a collector and those 
associated with the control building. 


Collector auxiliary device power requirements have been found Co be 

primarily a function of collector engine cooling requirements and for 

a collector in the 40-80 kW size range would he approximately 1.6 Id I 

• 

(per collector). 

Control building auxiliary device power requirements have been found to 
be a complex function of site environment, station operating schedule 
and energy storage requirements. Power requirements for these loads for 
the no storage, 643 kWh* storage (typical for 0.4 capacity factor station), 
7,S4S kWh* storage (typical for 0.7 capacity factor station) would be 
approximately (14 kW)/n, (32 kW)/n and 254 kW)/n, respectively, 
where "a” is the number of collectors in the collection field. 


Rev . 1 : 

1. Heat dissipation of building lights test equipment and outlets 
deleted from air conditioning equipment energy demand calculations 
based on plans for un-atiended station operation. 

2. Revised air conditioning equipment power requirements incorporated. 


* Due to system design changes, these figures for the P-75 engine are 
now 2,132 and 13,625 kWh, respectively. The storage values for the 
P-40 engine are 2,534 and 14,176 kWh, respectively. 


* 
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1.0 GENERAL CONSIDERATIONS 

Station auxiliaries, such as collector engine cooling fans, the power 
control microprocessor, tracking drive units and control building 
air conditioning equipment are necessary for station operation, but 
consume energy which must be subtracted from the gross station generating 
capacity to determine the net station output. These auxiliaries fall in 
two classes; i.e., those directly associated with the collectors and 
those associated with the control building and its installed equipment. 

Power requirements for the auxiliaries associated with any single collector 
are primarily a function of heat engine size and efficiency. Energy require- 
ments are a function of the foregoing factors and collector operating time, 
which has been assumed to average 10.14 hours per day, based on data 
provided in Table 1 of Reference l! Both energy and power requirements 
may be calculated in a fairly direct manner. 

Power and energy requirements for control building ancilliaries are 
recursive in nature and require assumption of a station operating schedule 
to calculate absolute values. ?or purposes of this report it has again 
been assumed that the period of direct energy collection will average 
10.14 hours per day. In those cases where storage is provided, it has 
also been assumed that the batteries will be recharged in the eight-hour 
period centered on the solar meridian and that they will be discharged 
through the inverter to the utility grid in the period immediately 
following the period of direct collection. Reasonable variations in the 
above described station operating schedule will not significantly affect 
the results presented. 


2.0 COLLECTOR AUXILIARIES 

Connected loads directly associated with a collector of 40-30 kW output 



rating are estimated to 

be as 

follows: 


p fp 

0.99 

kW 

(1 HP engine cooling fan & pump) 


Pec 

0.10 

kW 

(engine control system) 


Pbc 

0.05 

kW 

(engine starting battery charger) 


Pdc 

1.10 

kW 

(collector declination drive) 


Pha 

4.40 

kW 

(collector hour-angle drive) 


Ppc 

0.10 

lew 

(power control assembly) 

c 

? C 

6.74 

kW 

(collector total) 


i> 



* The more recent Barstow-based analysis shows that average operating 
time is approximately 8 hrs/day. 
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COLLECTOR AUXILIARIES continued: 

Energy requirements of each auxiliary device are a function of its power 
requirement and its operating duty cycle. The engine cooling fan & pump 
motor will ( for example) operate continuously while the collector is in 
service and for a 0.25 hour cool-down period iraiediately following the 
period of operation. This would give a duty cycle of (10.14 + 0.25)/24 
or 0.4329. Duty cycles for other auxiliary devices have been calculated 
on the basis of estimated device utilization and tabulated below with 
the rcsulatant energy requirements 





Power (TcW) 


Durv Cvcle 


Ers/Tav 

Snersrv 

*«? 

9 

0.99 

X 

0.4329 

X 

24 

10.29 

Zee 

9 

0.10 

X 

0.4225 

X 

24 

1.01 

Zbc 

9 

0.05 

X 

1.0000 

X 

24 

1.24 

Zdc 

9 

1.10 

X 

0.00758 

X 

24 

0.20 

Zha 

9 

4.40 

X 

0.0104 

X 

24 

1.10 

Spc 

9 

0.10 

X 

0.4225 

X 

24 

1.01 

Total 

Average Daily Energy Per Collector 


Zc 

14.35 


For an average dally collection period of 10.14 hours and a generator 
efficiency of 0.9, the average power required from the output shaft of 
each collector engine would be: 

? ' a I ^-'rn ^ m l 63 

(10. 14) (0.90) L '°- > 


Since the energy requirement of the major user (Efp) is nearly directly 
proportional to collector operating time, P c ' will be nearly independent 
of operating time and thus will not vary significantly throughout the year 
For smaller collector engines with an output rating in the range of 20-40 
ktf, all auxiliary loads would be the same except the engine cooling fan 
and pump motor which would be proportionally sized down to 1/2 E? and 
would draw 0.59 ktf. 3y calculations similar to the above, it can be shown 
that the average power requirement would be 1.17 kW per collector. 


COLLECTOR AUXILIARIES continued 


For larger collector engines with an output rating in the range of 30-120 kW, 
all auxiliary loads would again be the same except the engine cooling fan 
and pump motor which would be proportionally sized up to 2 HP and would 
draw 2.00 kW. Again, by calculations similar to above, it can be shown 
that the average power requirement would be 2.78 kW per collector. 

3.0 CONTROL 3LTLDETG AUXILIARIES 

Power requirements for control building auxiliary loads are a strong function 
of energy storage requirements, which in turn are a function of required 
station capacity factor and the percentage of capacity factor which can 
be satisfied from direct generation. It is therefore necessary to consider 
three major cases; i.e., the no storage case, the 0.4 capacity factor case 
and the 0.7 capacity factor case. 

For all three cases, it has been assumed that a well-constructed and well- 
insulated building, having 60#/sq. ft. wall construction, 20#/sq. ft. roof 
construction and an insulation "U" factor of 0.1, will be provided to house 
the station power control, storage (where used) and distribution equipment. 
Design ambient conditions have been assumed to be as follows, based on 
technical equipment requirements and typical sice conditions, as delineated 
in the Barstcv weather data. 


Inside Temperature: 

25° ± 

8°C 

Outside Summer : 

46 . 1°C 

dry bulb 


25.0OC 

wet bulb 

Outside Vinter : 

1.7°C 

dry bulb 

Daily Range : 

20.7°C 



The general calculation procedure used has been to assume that (for each case) 
a building of the configuration shown in the 3rd Progress Review Meeting 
will be oriented on the site to minimize solar exposure. Solar and trans- 
mission heat gain has then been calculated and totalled for the four wall 
and roof exposures (it has also been averaged over the floor area as a 
convience for rapidly considering the effect of changes la building size). 


CCNTHCL 3LTLDCTG AUXILIARIES continued: 


Heat dissipation has then been calculated and tabulated for all station 
auxiliaries in each room. Heat dissipation from the main power switchboard 
and storage subsystem components (where applicable) for the various modes 
of operation have then been calculated, tabulated and sussed with the 
aforementioned heat loads to compute the required sire ( 3 ) for the air 
conditioning equipment. Air conditioning equipment power requirements are 
then determined for each mode of operation based on the percentage of the 
installed capacity actually being used. This power is added to the other 
power flowing through the main switchboard and heat dissipation from the 
switchboard is re-calculated and re^totalled with the other heat loads 
to determine an adjusted value for building heat load and percentage of 
air conditioning capacity in use and hence a duty cycle for the air 
conditioning equipment. Although further Iterations would improve the 
mathematical precision of the result, the single iteration undertaken 
assures that the accuracy of the result will be determined by the accuracy 
of the input data assumptions. 

3.1 No Storage Case: 

A building of approximately 630 square feet (20' x 34') will be required 
to house the power control & distribution equipment, a small maintenance 
3hcp and a rest room. 3ased on the above stated design conditions, build- 
ing heat loads would be as follows: 


Bldg. Heat Gain 

680 sq. ft. 

X 

13.0 

3tuh/ sq.ft. 

m 

10,200 Btuh 

Lighting 

1, 300W 

X 

3.41 

Btu/watt hr. 


4,433 

* 

Outlets 

1,200W 

X 

3.41 

3tu/watt hr. 

m 

4,092 

★ 

Control equip. 

1, 900W 

X 

3.41 

Btu/watt hr. 

m 

6,479 


Main Switchboard 

3, 13 1W 

X 

3.41 

Btu/watt hr 

SB 

10, 677 


Peop le 

2 

X 

215 

Btuh/person- 

m 

430 

* 


Total Sensible Heat 36,311 Stub 

* Used for air conditioning unit capacity calculation, but not for energy 
calculation due unattended nature of station operation. 
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No Storage Case continued: 


o 




The control room cooling requirement could (typically) be satisfied by 
a Pomona Air ACC-300 cooling unit and a 3ohn RDD-5 air-cooled condenser. 
The cooling requirement in each of the smaller rooms could be satisfied 
by a Carrier iCC0051. The total connected load of all of the above units 
would be 14.10 kW. 

Connected loads in the control building are estimated to be as follows: 


?np 

m 

0.76 

kW 

(microprocessor) 

Pet 

9 

0.14 

kW 

(CRT terminal) 

Pad 

m 

1.00 

kW 

(A-D converters) 

Pel 

9 

0.30 

kW 

(control room lights) 

Pnii 

m 

0.40 

k W 

(maintenance shop lights) 

Prl 

m 

0.10 

k W 

(rest room lights) 

Pto 

m 

1.20 

kW 

(test equip. & outlets) 

Pac 

9 

14.10 kW 

(air-conditioning equip.) 

*3 

9 

18.50 


(building total) 


Znergy requirements for the above loads are a function of the power require- 
ment of each load and its operating duty cycle. In this case, it has been 
assumed that the microprocessor will operate 24 hours per day and chat the 
other technical loads will operate 10.14 hours out of the 24-hcur day. 

Lighting is assumed to be required an average of 4 hours per week and 
test equipment and outlet power is assumed to be required 2 hours per week. 

The air conditioning equipment will operate. 567. of the time during 10.14 
operating day and will be shut down during the nonoperating period to conserve 
energy. Under these conditions energy requirements would be as follows: 




Power (V.O 


Du tv Cvcle 


Hrs/Dav 


Enersv Rents. 


9 

0.76 

X 

1.0000 

X 

24 

9 

13.24 kWh 

Set 

9 

0.14 

X 

0.4225 

X 

24 

9 

1.42 

Sad 

9 

1.00 

X 

0.4225 

X 

24 

9 

10.14 

Eel 

9 

0.80 

X 

0.0238 

X 

24 

9 

0.46 

Eml 

9 

0.40 

X 

0.0238 

X 

24 

9 

0.23 

Erl 

9 

0.10 

X 

0.0238 

X 

24 

9 

0.06 

Eto 

9 

1.20 

X 

0.0119 

X 

24 

9 

0.34 

Sac 

9 

14.10 

X 

0.2789 

X 

24 

9 

94.38 



Total Daily 

Energy Requir 

ement 

-3 

9 

125.27 kwh 
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So Storage Case continued: 


The average power requirement for Che building auxiliary equipment loads 
during a 10.14 hour operating day would be: 



125.27 

10.14 


12.35 


If there are ’a' collectors in the collection field, the average auxiliary equipment 
power required at the engine output shaft of each collector during Che operating 
day would be: 


V 


12 • 35 

a (0.997) (0.981) (0.900) 


14.03 


a 


xK 


The above requirement is additive Co the requirement for collector auxiliary 
equipment power requirements defined in paragraph 2.0. It includes no 
allowance for perimeter (or other) site security lighting. 


2 


0.4 Capacity Factor Storage Case: 

Power requirements for control building auxiliary equipment loads for Che 
0.4 capacity factor case are a strong function of Che stored energy require- 
ment, which is equal to the total energy requirement less the amount of 
energy available from direct generation. Computations in this section assume 
that a capacity factor of 0.373 will be achieved from direct generation 
per Table 4 of Reference 1. Any significant variations in this value 
will affect required battery ampere-hour capacity, building size, collector 
and inverter operating time, air conditioning equipment size and energy 
requirements. They will have a significant effect on building related 
auxiliary equipment power and energy requirements. Average daily stored 
energy delivery to the utility grid for a station having a direct generation 
capacity factor of 0.373 would be: 


2S 


tl ; 0 K10V .37S0 ) (0.400-0. 373) . M ^ 

365 


The total storage capacity required to deliver this amount of energy would be: 


E batt 


643 

HV) C)T) (77 L) 


64a • 

(0.999) (0 . 98) (0 . 94 ) 


704.1 kWh 


wtaere^ ^ j and^f are the efficiencies of the HV connection equipment, 
LV-EV transformer and inverter, respectively. 


▲ building to house the storage subsystem of this size, in addition to the 
basic power control and distribution equipment, a maintenance shop and rest 
room would require a floor area of approximately 3,072 square feet. 
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3.2 0.4 Capacity Factor Storage Case continued: 


( 


T\ 


Technical equipment heat loads within the building will vary significantly 
during a typical day's operation depending on whether the storage subsystem 
is in a float, charge, or discharge node. 

In the float charge node, which last 2.14 hours of the typical operating 
day, the building heat loads would be as follows: 


In the charge mode, which would last 3.0 hours of the typical operating 
day., the basic building heat loads would be increased by the converter 
and battery heat losses. The total sensible heat load for this mode of 
operation is summarized below. 


Bldg . Heat Gain 

3,072 sq. f 

t. X 

10.0 

Btufa/sq.ft. 

* 

30,720 

3 tun 

Lighting 

3, 900W 

X 

3.41 

3tu/watt hr. 

3 

13,299 

* 

Outlets 

1.50CW 

X 

3.41 

Btu/watt hr. 

3 

5,115 

* 

Control Equipment 

1.900W 

X 

3.41 

Btu/watt. hr. 

a 

5,479 


Main Switchboard 

3,414W 

X 

3.41 

Btu/watt hr. 

a 

11,642 

• 

Converters 

3.200W 

X 

3.41 

Btu/watt hr. 

a 

10, 912 


Inverters 

40 CW 

X 

3.41 

Btu/watt hr. 

a 

1,364 


Battery 

10.360W • 

X 

3.41 

Btu/watt hr. 

a 

35,328 


Pe' >le 

2 

X 

215 

Btuh/person 

a 

430 

* 


• Total 

Sensible 

Heat 



115,239 

3 tun 


* Used for air conditioning unit capacity calculation, but not for energy 
calculation due unattended nature of station operation. 
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Bldg. Heat Gain 

3,072 sq.ft 

X 

10.00 Btuh/ sq.ft. 

a 

30,720 Btuh 

Lighting 

3, 900W 

X 

3.41 Btu/watt hr. 

m 

13,299 * 

Outlets 

1, 500W 

X 

3.41 3tu/vatt hr. 

a 

5,115 * 

Control Equipment 

1, 900W 

X 

3.41 Btu/watt hr. 

a 

6,479 

Main Switchboard 

3,071W 

X 

3.41 Btu/watt hr. 

a 

10,472 

Converters 

40 OW 

X 

3.41 Btu/watt hr. 

a 

1,364 

Inverters 

400W 

X 

3.41 3tu.watt hr. 

a 

1,364 • 

Battery 

0 

X 

3.41 Btu/watt hr. 

a 

0 

People 

2 

X 

215 Btuh/person 

a 

430 * 


Total 

Sensible Heat 


69,243 3tuh 








A — 
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3.2 0.4 Capacity Factor Storage Case c ontinued: 

Is Che discharge mode, which would lasc 0.93 hours of a cynical operating 
day, the inverter and storage battery losses would contribute to the 
total sensible heat. Air conditioning requirements during this period 
would therefore be as follows: 


31dg. Heat Gain 

3,072 sq. ft. 

X 

10.0 

3 tun/ sq.ft. 

m 

30,720 Btufa 

Lighting 

3,900V 

X 

3.41 

3tu/vatt hr. 

m 

13,299 * 

Cutlets 

1,500W 

X 

3.41 

3tu/vatt hr. 

m 

5,115 * 

Control Equipment 

I, 900V 

X 

3.41 

3tu/vatt hr. 

m 

6,479 

Main Switchboard 

2,422V 

X 

3.41 

Btu/watt hr. 

m 

8,259 

Converters 

400V 

X 

3.41 

3tu/watt hr. 

m 

1,364 

Inverters 

45,640V 

X 

3.41 

Btu/watt hr. 

9 

155,632 

Battery 

45, 120V 

X 

3.41 

Btu/watt hr. 

9 

153,859 

People 

45, 120V 

■ X 

215' 

3tuh/person 

- 

430 * 


Total Sensible Heat 375,157 3tuh 

As can be seen from the above tabulations, there is an approximate 5:1 
variation in total heat load encompassed by the various operating modes. 
Both the control room and battery room cooling requirements could be 
satisfied by separate air conditioning systems consisting of a Pomona 
ACC-500/3ohn RDD-9 combination supplemented by two Trane SACA 754 units, 
all operating under staged thermostatic control. Cooling requirements of 
the two smaller rooms could be satisfied, as before, by two Carrier KX0051 
units. 3ased on the foregoing, connected loads in the control building are 
estimated to be: 


?np 

m 

0.76 kV 

(microprocessor) . 

?ct 

m 

0.14 

(CRT terminal) 

Pad 

m 

1.00 

(A-D converters) 

?cl 

m 

1.30 

(control room lights) 

Pbl 

m 

1.60 

(battery room lights) 

Pml 

m 

0.40 

(maintenance shop lights) 

Prl 

m 

0.10 

(rest room lights) 

Pto 

m 

1.50 

(test equip. & outlets) 

Pac 

m 

93.09 

(air conditioning equip) 

?B 

9 

100.39 kW 

(building total) 
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0.4 Capacity Factor Storage C3se continued; 

2nergy requirements for the above loads are, again, a function of 
the pover requirement and duty cycle of each load. Based on an assumed 
mode of station operation which would result in direct generation and 
power delivery to the grid for 10.14 hours, battery charging for 3.00 
hours and delivery of power from storage to the grid for 0.93 hours, 
the average daily auxiliary equipment energy requirement would be as 
follows: 




Rower 

OcW) 

Du tv Cvcle 


Hrs/Dav 

Znerzv Rants 

2^ 

a 

0.76 

X 

1.0000 

X 

24 

- 18.24 kWh 

Set 

a 

0.14 

X 

0.4225 

z 

24 

- 1.42 

Sad 

a 

1.00 

X 

0.4225 

X 

24 

- 10.14 

Bel 

a 

1.80 

X 

0.0238 

X 

24 

- 1.03 

Sbl 

a 

1.60 

X 

0.0238 

X 

24 

- 0.91 

Sal 

m 

0.40 

X 

0.0238 

X 

24 

- 0.23 

Sri 

a 

0.10 

X 

0.0238 

X 

24 

- 0.06 

Eto 

a 

1.50 

X 

0.0119 

X 

24 

- 0.43 

Sacf 

a 

93.09 

X 

0.0102 

X 

24 

• 22.79 

Saco 

a 

93.09 

X 

0.0731 

X 

24 

. 163.32 

®acd 

a 

93.09 

X 

0.0314 

X 

24 

- 70.15 


Total Daily Energy Requirement 


288.72 kWh 


The average auxiliary equipment power requirement for the building loads 
during a 10.14 hour operating day would be: 


V 


238.72 

10.14 


28.47 kW 


If there are 'V collectors in the collection field, the average power 
required at the engine output shaft of each collector would be: 


pb" 

The above 
equipment 
allowance 


28.47 32.34 ... 

n(0.997) (0.981) (0.90) “ n 

requirement is additive to the requirement for collector 

power requirements defined in paragraph 2.0. It includes 

for perimeter (or other) site security lighting. 


auxiliary 

no 
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0.7 Capacity Factor Storage Case; 

Power requirements for control building auxiliary equipment loads for the 
0.7 capacity factor case are also a strong function of the stored energy 
requirement, which is equal to the total energy requirement less the amount 
o£ energy available from direct generation. 


Computations in this section assume that a capacity factor of 0.373 will 
be achieved from direct generation per Table 4 of Reference 1. Any signi- 
ficant variations in this value will affect required battery ampere-hour 
capacity, building site, collector and inverter operating time, air condition- 
ing equipment size and energy requirements . They will have a significant 
effect on building related auxiliary equipment power and energy requirements. 
Average daily stored energy delivery to the utility grid for a station 
having a direct generation capacity factor of 0.373 would be: 

Is . (1.0MX03K3760K0.700-0.373) . ^ ^ 

The total storage capacity required to deliver this amount of energy 
would be: 

7.34 3 7.343 . 

Ebatt - HV) y) ( 7 ) x) " (0.999)(0.98)(0.94) " 8 ’ 52S 


A building to house the storage subsystem of thi3 size, in addition to 
the basic power control and distribution equipment, a maintenance shop 
and rest room would require a floor area of approximately 3,443 square 
feet. Technical equipment heat loads within the building will vary 
significantly during a typical day's operation depending on whether 
the storage subsystem Is in a float, charge, or discharge mode. 


In the float charge mode, which would last 2.14 hours of the typical 
operating day, the building heat loads would be as follows: 




k 

r 


• • 

3.3 

0.7 Caoacitv Factor 

Storaee Case 

continued 




1 

§ <1 ) 
6 - 

Bldg. Heat Gain 

8,448 sq.ft. 

X 

3.50 

3 tun/ sq.ft. 

71,808 

3tuh 

Lighting 

8,500V 

X 

3.41 

Btu/watt hr. 

28,985 

★ 


Outlets 

1, 3C0W 

X 

3.41 

3tu/wact hr. 

6,138 

it 


Control Equipment 

1,900V 

X 

3.41 

3tu/watt hr. 

• 6,479 



Main Switchboard 

3, 152V 

X 

3.41 

3tu/watt hr. 

10,748 



Converters 

400V 

X 

3.41 

Stu/watt hr. 

1,364 



Inverters 

400V 

X 

3.41 

Btu/watt hr. 

1,364 


• 

3attery 

0 

X 

3.41 

Btu/watt hr. 

0 



People 

2 

X 

215 

3tuh/person 

430 

it 



■( 


r 


Total Sensible Heat 


127,316 3tuh 


la the charge node, which would last 8.0 hours o£ the typical operating 
day, the basic building heat loads would be increased by the converter 
and battery heat losses. The total sensible heat load for this node 
of operation is submarined below. 


Bldg Heat Gaia 

8,448 sq.ft. 

X 

8.5 

Btuh/ sq.ft. 

m 

71,808 

3tuh 

Lighting 

8,500V 

X 

3.41 

Btu/watt hr. 

9 

23,985 

* 

(Xitlets 

1,300V 

X 

3.41 

Btu/watt hr. 

m 

6, 138 

★ 

Control Equipment 

1,900V 

X 

3.41 

Btu/watt hr. 

9 

6,479 


Main Switchboard 

7,513V 

X 

3.41 

Btu/watt hr. 

9 

25,619 


Converters 

38,730V 

X 

3.41 

3tu/watt hr. 

9 

132,240 


Inverters 

40CV 

X 

3.41 

Btu/watt hr. 

9 

1,364 


Battery 

125.40CV 

X 

3.41 

3tu/watt hr. 

9 

427,614 


People 

2 

X 

215 

Btuh/person 

9 

430 

* 


Total 

Sensible 

Heat 


700,677 

Btuh 

In the discharge mode, which would last 11 

.21 hours of a 

typical operating 


day, the inverter and storage battery losses would contribute to the 
total sensible heat. Air conditioning requirements during this period 
would therefore be as follows: 

* Used for air conditioning unit capacity calculation, but not for energy 
calculation due unattended nature of station operation. 
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0.7 Capacity Factor Storage Case continued: 


31dg Heat Gain 

8,448 sq.ft. 

X 

3.5. 

3tuh/sq.ft. 

a 

71,308 

3 tun 

Lighting 

3,5009 

X. 

3.41 

Btu/vatt hr. 

a 

28,985 

* 

Outlets 

1 , sccw 

X 

3.41 

3tu/watt hr. 

a 

6,138 

* 

Control Equipment 

1,9009 

X 

3.41 

3tu/watt hr. 

M 

6,479 


Main Switchboard 

2,5139 

X 

3.41 

3tu/vatt hr. 

a 

3,569 


Converters 

4009 

X 

3.41 

3tu/watt hr. 

SB 

1,364 


Inverters 

45,6409 

X 

3.41 

Btu/vatt hr. 

a 

155,632 


Battery 

45, 1209 

X 

3.41 

3tu/vatt hr. 

a 

153,359 


People 

2 

X 

215 

3tuh/person 

M 

430 

* 


Total 

Sensible 

Heat 


433,264 

3 tun 


As can be seen from Che above tabulations, there is an approximate 6:1 
variation in total heat load encompassed by the various operating modes. 
The control room cooling requirement could be satisfied by two Pomona 
ACC -800-2 /3ohn RDD-L6 combinations. The battery room would require three 
Pomona ACC-L500-2/3ohn SDD-27 combinations. All units would operate 
under staged thermostatic control. Cooling requirements of the two 
smaller rooms could be satisfied, as before, by two Carrier ££0051 units. 
3ased on the foregoing, connected loads in the control building are 
estimated to be: 


?mo 

SB 

0.76 kW 

(microprocessor) 

?ct 

a 

0.14 

(CUT terminal) 

Pad 

SB 

1.00 

(A-D converters) 

»cl 

a 

1.80 

(control room lights) 

Pbl 

a 

7.20 

(battery room lights) 

P®1 

m 

0.40 

(maintenance shop lights) 

Prl 

a 

0.10 

(rest room lights) 

Pto 

a 

1.80 

(test equip. & outlets) 

Pac 

a 

175.73 

(air conditioning equip) 

P B 

a 

188.93 

(building total) 


* Used for air conditioning unit capacity calculation, but not for 
energy calculation due unattended nature of station operation. 
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icitv Factor Storage Case continued: 


Energy requirements for the above loads are, again, a function of 
the power requirement and duty cycle of each load. 3ased on an assumed 
mode of station operation which would result in direct generation and 
power delivery to the grid for 10.14 hours, battery charging for 3. CO 
hours and delivery of power from storage to the grid for 11.21 hours, 
the average daily auxiliary equipment energy requirement would be as 
follows: 




Power (W) 


Durr Cvcle 


Hrs/Dav 


Enerzv Rents 

Srrrp 

9 

0.76 

X 

1 0000 

X 

24 

m 

18.24 kWh 

Set 

m 

0.14 

X 

0.4225 

X 

24 

9 

1.42 

Sad 

m 

1.00 

X 

0.4225 

X 

24 

m 

. 10.14 

=01 

m 

1.30 

X 

0.0238 

X 

24 

m 

1.03 

2bl 

9 

7.2 

X 

0.0238 

X 

24 

m 

4.11 

Eml 

u 

0.40 

X 

0.0238 

X 

24 

m 

0.23 

2rl 

9 

0.10 

X 

0.0238 

X 

24 

9 

0.06 

2to 

m 

1.80 

X 

0.0119 

X 

24 

9 

0.51 

3 acf 

m 

175.73 

X 

0.0113 

X 

24 

9 

47.66 . 

B ace 

m 

175.73 

X 

0.2624. . 

X 

24 

m 

1106.67 

2acd 

m 

175.73 

X 

0.2565 

X 

24 

- 

1081.79 



Total Daily Energy Requirement Eg' 

9 

2271.86 kWh 


The average auxiliary equipment power requirement for the building loads 

during a 10.14 hour operating day would be: 

2271.86 

n_» . a 224.05 

10.14 

If there are "a" collectors in the collection field, the average power 
required at the engine output shaft of each collector would be: 

224.05 _ 254.53 

53 * n(0. 997) (0.981) (0.90) " a ‘ 

The above requirement is additive to the requirement for collector auxiliary 
equipment power requirements defined in paragraph 2.0. It includes no 
allowance for perimeter (or other) site security lighting. 
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3.4 


Energy Requirements for Security Lighting; 

Perimeter security lighting (if required) for the 0.4 capacity factor 
station would typically consist of 32 175W halogen lamps for which 

power requirements would be: 


PL * 32 X 0.175 - 5.6 kW 

And for an annual average daily sunset to sunrise period of 11. S3 hours 
would require: 

El * 5.6 X 11.93 - 66.3 kWh 


If there "a 1 * collectors in the collection field, the average power 
required at the engine output shaft would be: 

6fe.3 7.48 


V 


n(0 . 997) (0 . 981) (0.90) (10 . 14) 


kU 


3.5 Cold Weather Heating Requirements: 

Cold weather heating has not been specificially addressed in the preceding 
paragraphs, since excess aeac will be available during all periods of 
equipment operation and all equipment items are suitable for a non- 
operating temperature environment of 0°C, which is less than the minimum 
listed in the weather data for the typical site location. 

Heat balance calculations should, however, be made for sites having a 
lower temperature irtnimnn to determine whether heat dissipation from the 
installed equipment and re-radiation from the thermal mass of the building 
structure will maintain the building interior above 0°C or whether 
supplemental heating will be required. 


REFERENCES 

1. Pons, R.L., ’’Preliminary System Performance Analysis”, Technical 
Report No. SPS-007, dated 27 September 1978. 
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LIFE CYCLE COST EQUATIONS* 


This Appendix presents a summary of the life cycle cost (LCC) equations used 
for the system energy cost analyses. An effective energy "price" is required 
in order to recover the costs associated with the purchase, installation, 

operation and maintenance of a utility-owned solar electric system. This 

energy "price" is expressed in t erms of the levelized busbar energy cost (BBEC) 
The methodology for determining BBEC is presented below. 


(a) Present Value of Capital Investment (Cl ) 

. n * «.? 


(1 + r t TT 


where: 


CI c = Capital outlay in year t 

g c * Escalation rate for capital costs (0.06) 

k * Cost of capital (0.086) 

p - Y - Y 

r co p 

Y cq * First year of commercial operation (1988 for SPS example) 

Yp ■ Price year for cost information (1978) 

j = Y - Y +1 
J t co 

Y c * Year of given investment outlay 

(b) Annualized System Resultant Cost, AC 

(without fuel costs or operating costs) 

AC = (l+g)~ d [FCR • Cl + CRF, M MNT ] 

pv k,N pv 

* 0.08739 Cl +0.04773 MNT 

pv pv 


*0btained from: "The Cost of Energy from Utility Owned Solar Electric 

System", JPL Report 5040-79, ERDA/JPL 1012-76/3, June 1976 


H-l 


where: 


g « Rate of general inflation (0.06) 

d « Y - Y. - 1988 - 1978 =- 10 

co b 

Y, = Base vear = 1978 

D 

FCR * Annualized Fixed Charge Race (0.1565) 

CRF. .. = Capital Recovery Factor (0.0939) 

MNT ■ Present value of recurring maintenance costs 
pv 


(c) Present Value of Maintenance Costs, MNT pv for uniform growth. 


MNT pv =■ (1 + 


_1 + 8m' 

: 

!'■ 

" l + 

m 

s ) 

i 

B 

t>0 

1 

1 


l ♦ :< 

j 


X 


47.274 X 


where: 


g ^ » Escalation race for maintenance costs (0.07) 

N ■ System operating lifeline (30 years) 

X = Recurrent costs ($) 
o 

for non-uniform growth. 


MNT 




1 + 8 , 


-|J 


m 

1 + K 


X 


1.96715 (0.98527) J X 


(3) 


(4) 


(d) Levelized Busbar Energy Costs, BBEC 


BBEC 


AC 

MWH 

A 


(5) 


H-2 


where: 


MWH * ACF x 8760* = actual plant output 

A 

ACF = Annualized capacity factor, including allowance for 
unscheduled maintenance 


BBEC_ * 9.9759(10) ° AC (for leap vear) 
C c 


0M * 5.4491(10) • AC (for leap year) 


Due to capital expenditures 
Due to maintenance expenditures 


BBEC = 9. 9759 (10) ' 6 AC + 5 . 4491(10) “ 6 AC„ 

C M 


* ■ 8784 for leap year 
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T.R. No. SPS-20 
(Revision A) 


RELIABILITY PREDICTIONS FOR THE 
3ASELINE SYSTEM AND OTHER SYSTEMS 


1 . INTRODUCTION 

This cechnical report (TR) revises the predicted cumulative mean- 
time-betveen-failure (CMTBF), cumulative mean-time-to-repair (CMTTR) and 
availability of full rated output during *-he scheduled power delivery 
time for the baseline Solar Small Power System (SPS) having a plant size 
of 1.0 megawatt electrical (MWe) output with an annual capacity factor 
(ACF) of 0.4. 

This TR also includes the predicted CMT3F, CMTTR and availability 
for additional plant sizes of 0.5 and 10.0 MWe with ACF 1 s of 0.36, 0.4 
and 0.7 comprising eight additional Solar SPS's. 

2. PREDICTION SUMMARY 

Availability is defined as the probability that a system will operate 
as required by its design at any point in time and considers only operat- 
ing time (uptime) and downtime, thus excluding idle time. Downtime refers 
to the corrective (unscheduled) maintenance action time required to return 
the system, subsystem, equipment or item to the operational state. Availa- 
bility, derived by the ratio of the system CMTBF to the sum of system 
CMTBF and CMTTR, corresponds to the uptime probability. 

The data presented herein are structured primarily to focus on the 
areas where corrective (unscheduled) maintenance will be required during 
the solar SPS operation. Table I summarizes the results of the reliability 
and availability prediction for the baseline system as well as the eight 
other systems considered. 

The data presented in Table II summarize the predicted subsystem 
failures and the corresponding subsystem failure distribution in percent 
to identify the subsystems where failures are most likely to occur over 
the system power output duration, T q . 

The data presented in Table III summarize the predicted cumulative 
subsystem downtime in hours resulting from system operation over the time 
duration, T q . The downtime distribution in percent for each subsystem is 
presented to identify the subsystem that drastically affects system 
availability or uptime. 

3 . CONCLUSIONS 

The following conclusions are based on data presented in Tables I, 

II and III. 
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I 


3. 


CONCLUSIONS - (concinued) 


• The cumulative system failure and the system downtime 
increase with the solar SPS complexity (ref. Table I). 

• The system complexity has definite effect on the CMT3F 
(ref. Table I). The CMT3F for the baseline system is 
222 hours. 

• The availability is relatively insensitive to the ACF 1 s 
considered for a given plant size (ref. Table I). The 
availability of the baseline system is 0.9868. 

• The collector and power subsystems of systems with ACF 
of 0.36 have the highest failure percentages due to the 
absence of power storage capability (ref. Table I). 

• The collector and power conversion subsystems of any 
solar SPS will contribute to more than 60 percent of the 
annual failures (ref. Table II). For the baseline system, 
the collector and power conversion subsystems will fail 
approximately 10 times during the year, corresponding to 

70.3 percent of the annual number of expected system failures. 

• More than 85 percent of the system downtime is due to the 
collector and power conversion subsystems (ref. Table III). 

For the baseline system, its collector and power conversion 
subsystem dowtime is approximately 38 hours per year, 
corresponding to 89.9 percent of the annual system downtime. 

• The prediction data indicate that ease of maintenance, 
repair or replacement must be designed into the system to 
reduce downtime. Modularization of items with relatively 
high failure rates and downtime should be considered to 
enhance accessibility, handling and quick replacement in 
the event of a failure. 

• Provision to switch-in a standby subsystem should be a con- 
sideration in the design to negate the effects of long down- 
time due to repair or replacement. The down t 1ms Zui shis con- 
sideration is the time it takes the subsystem to switch-in 

and to attain operational status. This approach can 
drastically increase subsystem availability in areas where 
the frequency of failure occurrence and the corresponding 
downtime are relatively high. 

4. ANALYSIS 

The solar SPS analyses are based on information provided by Enginee 
ing on nine different system configurations encompassing plant sizes of 
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--SYSTEM 

Plant 

Size 

DESCRIPTION— 
ACF Storage 

TABLE I . SYSTEM 

Power 

Output 

Duration 

(Hours) 

RELIABILITY AND 

rt VAI LABILITY 

PREDICTION SUMMARY 



Cumulative 

Failures 

Downtime 

(Hours) 

’ " “ (KCiUit 1 1UN- - - - - - 

CMTBF 

(Hours) 

CMTTR 

(Hours) 

Full Rated 
Avai lability 

0 . 5MWe 

0.36 

No 

2838 

6.09 

19.51 

466 

3.20 

0.9931 

0. 5MWe 

0.40 

Yes 

3168 

8.03 

22.85 

395 

2.85 

0.9928 

0. 5MWe 

0.70 

Yes 

5544 

17.61 

47.82 

315 

2.72 

0.9914 

l.OMWe 

0.36 

No 

2838 

11.06 

37.63 

257 

3.40 

0.9869 

* 1 • OMWe 

0.40 

Yes 

3168 

14.29 

42.32 

222 

2.96 

0.9868 

l.OMWe 

0.70 

Yes 

5544 

32.53 

91.54 

170 

2.81 

0.9837 

10. OMWe 

0.36 

No 

2838 

100.8 

366.6 

28 

3.64 

0.8849 

10. OMWe 

0.40 

Yes 

3168 

131.9 

412.5 

24 

3.13 

0.8846 

10. OMWe 

0.70 

Yes 

5544 

305.9 

896.4 

18 

2.93 

0.8600 


* Baseline System Calculations: 


CMTBF 

= 3168 Hours = 

14.29 Failures 

222 hours between failures 

CMTTR 

= 42.32 Hours 

14.29 Failures 

2.96 hours 

to repair 

Availability 

= CMTBF 

222 

= 0.9868 

CMTBF + CMTTR 

222 + 2.96 



TABI.E II. PREDICTED SUBSYSTEM FmtiURES AND FAILURE DISTRIBUTION 
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SUBSYSTEM FAILURES AND DISTRIBUTION 


Pl ant Power Transport/ Computer System 


Size 

ACF 

Data 

Col lector 

Conversion 

Distribution 

Control 

Storage 

Total 

0 . 5MWe 

0.36 

Failures 

2.32 

2.45 

0.27 

1.05 


6.09 



Distribution 

38.1% 

40.3% 

4.4% 

17.2% 


100% 

0. 5MWe 

0.40 

Failures 

2.56 

2.72 

0.36 

1.17 

1.22 

8.03 



Distribution 

31.8% 

33.9% 

4 . 5% 

14.6% 

15.2% 

100% 

0 . 5MWe 

0.70 

Failures 

5.39 

5.71 

1.14 

2.06 

3.31 

17.61 



Distribution 

30.6% 

32.4% 

6.5% 

11.7% 

18.8% 

100% 

l.OMWe 

0.36 

Failures 

4.62 

4.90 

0.49 

1.05 

— 

11.06 



Distribution 

41.8% 

44.3% 

4.4% 

9.5% 

“ “ “ 

100% 

* 1 . OMWe 

0.40 

Failures 

4.87 

5.17 

0.63 

1.17 

2.45 

14.29 



Distribution 

34.1% 

36.2% 

4.4% 

8.2 % 

17.1% 

100% 

l.OMWe 

0.70 

Failures 

10.52 

11.16 

2.19 

2.06 

6.60 

32.53 



Distribution 

35.3% 

37.5% 

4.2% 

0.9% 

22.1% 

100% 

10. OMWe 

0.36 

Failures 

46.2 

49.0 

4.6 

1.0 

— 

100.8 



Distribution 

45.8% 

48.6% 

4.6% 

1.0% 

* “ ™ 

100% 

10 . OMWe 

0.40 

Failures 

48.7 

51.7 

5.8 

1.2 

24.5 

131.9 



Distribution 

36.9% 

39.2% 

4.4% 

0.9% 

18.6% 

100% 

10. OMWe 

0.70 

Failures 

105.2 

111.6 

21.0 

2.1 

66.0 

305.9 



Distribution 

34.4% 

36.5% 

6.9% 

0.7 % 

21.5% 

100% 





*Baseline System 






TABLE III 


PREDICTED SUBSYSTEM DOWNTIME AND DISTRIBUTION 


SYSTEM SUBSYSTEM DOWNTIME (DT) AND DISTRIBUTION 

Plant Power Transport/ Computer System 


Size 

ACF 

Data 

Col lector 

Conversion 

Distribution 

Control 

Storaee 

Total 

0. 5MWe 

0.36 

DT (Hours) 

8.29 

9.80 

0.37 

1.05 

... 

19.51 



Distribution 

42.5% 

50.2% 

1.9% 

5.4% 


100% 

0 . 5MWe 

0.40 

DT (Hours) 

9.09 

10.88 

0.43 

1.17 

1.28 

22.85 



Distribution 

39.8% 

47.6% 

1.9% 

5.1% 

5.6% 

100% 

0. 5MWe 

0.70 

DT (Hours) 

19.16 

22.84 

0.93 

2.06 

2.83 

47.82 



Distribution 

40.1% 

47.8 % 

1.9% 

4.3% 

5.9% 

100% 

1 . OMWe 

0.36 

DT (Hours) 

16.50 

19.60 

0.48 

1.05 

— 

37.63 



Distribution 

43.8% 

52 . 1% 

1.3% 

2.8% 

“ “ ~ 

100% 

* 1 . OMWe 

0.40 

DT (Hours) 

17.35 

20.68 

0.56 

1.17 

2.56 

42.32 



Distribution 

40.1% 

48.9% 

1.3% 

2.8% 

6.0% 

100% 

1 . CMWe 

0.70 

DT (Hours) 

37.41 

44.64 

1.81 

2.06 

5.62 

91.54 



Distribution 

40.8% 

48.8% 

2.0% 

2.3% 

6.1% 

100% 

10. OMWe 

0.36 

DT (Hours) 

165.0 

196.0 

4.6 

1.0 

— 

366.6 



Distribution 

45.0% 

53.4% 

1.3% 

0.3% 

“ “ ~ 

100% 

10. OMWe 

0.40 

DT (Hours) 

173.5 

206.8 

5.4 

1.2 

25.6 

412.5 



Distribution 

42.1% 

50.1% 

1.3% 

0.3% 

6.2% 

100% 

10. OMWe 

0.70 

DT (Hours) 

374.1 

446.4 

17.6 

2.1 

56.2 

896.4 



Distribution 

41.7% 

49.8% 

2.0% 

0.2% 

6.3% 

100% 


* Baseline System 


4. ANALYSIS - (continued) 

0.5, 1.0 and 10.0 MWe, each plant size considered for ACF's of 0.36, 

0 . 4 and 0.7. 

4.1 RELIABILITY MODEL 

The basic reliability model of the solar SPS is configured as a 
series model illustrated by the block diagrams of Figure 1. The basic 
difference between Figure 1A and IB is the presence of the storage sub- 
system for the ACF's indicated. The main difference among the systems 
analyzed is the variation in complexity of the subsystem as summarized 
by Table IV for (N) number or quantity of item per subsystem. 

4.2 ASSUMPTIONS AND CONDITIONS 

The annual solar SPS operational or rated power output duration is 
determined by the ACF (see 4.3.1). Failure occurrence during this time 
duration will require corrective maintenance action to restore the 
system to full operational status by replacing, repairing or adjusting 
the component or subsystem which caused interruption of service. 

Preventive maintenance is performed at regular intervals when the 
SPS is not delivering scheduled power in order to maintain the system 
in a condition consistent with its designed level of performance, 
reliability, and where applicable, safety. 

Preventive maintenance typically involves servicing, inspections 
and minor or major overhauls during which 

a. regular care is provided to the normally operating subsystems 
and components which require such attention (lubrication, 
cleaning, adjustments, alignment, etc.), 

b. periodic functional check is performed to ascertain the 
integrity of the subsystem, 

c. marginal or erratic components and subsystems are repaired 
or removed and replaced, and 

d. components which are nearing wearout condition are replaced 
or overhauled. 

4.3 DEFINITIONS, EQUATIONS AND CALCULATIONS 

The following definitions and equations were established and used 
to calculate die various reliability parameters. 
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STORAGE 

SUBSYSTEM 


i 

oo 


A. Model applicable to the Solar SI’S with ACE of 0.4 and 0.7. 



COMPUTER 

CON'IKOL 

SUBSYSTEM 


B. Model applicable to the Solar SPS with ACE of 0.T6. 


Eigure T. Basic Reliability Model Block Diagrams. 


AX 












TABLE IV. PREDICTED FAILURE RATE SUMMARY 

BASED ON 

PLANT SIZE AND ACF COMPLEXITY 


4 a/m Aa/aju*)c dA^AC/r/ /zvcrex ('ac^) 


l.US 

*iu 


/f/#y i 


/*ns 


3 ms 


34H.S 




/Z* i 


KM C*\ 

















































4.3.1 ANNUAL CAPACITY FACTOR. 


The ACF is defined by che racio: 


ACF - Ja. 



where: T 0 * total time a given system delivers 

rated power output per year, and 

Tb * time base of 330 calendar days or, 

7920 hours per year. 

For ACF's of 0.36, 0.4 and 0.7, the corresponding T 0 's are 2838, 

3168 and 5544 hours, respectively. 

Or. a daily basis, the power delivery time breakdowns listed below 
relative to solar or storage source are assumed for the different ACF's. 
The time durations as shown form the basis of item operational time 
of Table V. 

Time (Hours)-- 


Solar Storage Total 


ACF 

Source 

Source 

Time 

0.36 

3.6 

None 

8.6 

0.40 

3.6 

1.0 

9.6 

0.70 

8.6 

8.2 

16.8 


4.3.2 FAILURE AND FAILURE RATE. Failure is defined as the inability 
of a part, a component or an equipment to perform as specified, 
or a system to continue operation as required, resulting in 
repair, removal or replacement action. 

Failure rate, the number of failure occurrences per unit time, 
is generally expressed in terms of failures per million operat- 
ing hours. The single item failure rates (X) tabulated in 
Table IV were estimated by using parts count, equipment simi- 
larity and complexity assessment techniques. The basic fail- 
ure data for the drive controls of the collector subsystem were 
provided by WDL Division of FACC, based on their field experiences. 
Other sources of the basic part/conrponent failure rates were the 
following: 


Part/Component Data Source 

a. Electronic MIL-HDBK-217B, "Military Standardization 

Handbook, Reliability Prediction of 
Electronic Equipment", dated 9-20-74; 
Notice 1, dated 9-7-76; Note 2, dated 
3-17-78. 
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4.3.2 FAILURE AND FAILURE RATE. - (continued) 


Part/Component 


Data Source 


b. Nonelectronic, 
mechanical, etc. 


NPRD-l, "Nonelectronic Parts Reliability 
Data", Summer 1978, ITT Research Institute/ 
Rome Air Development Center (RADC) 


c. Nonelectronic, 
mechanical, etc. 


RADC-TR-75-22, "Nonelectronic Reliability 
Notebook", dated January 1975, RADC. 


d. Mechanical Piece 
Part 


T- 70-48891-007, "Handbook of Piece Part 
Failure Rates", (HPPF) , dated 6-22-70, 
Martin Marietta, Denver Division, GIDEP 
Access No. E031-1273. 


The total item failure rate (X^) is defined by the expression 


where: N * quantity of a given item 

X = single item failure rate. 


The total expected item failure quantity per year is determined 
by the expression 


where: X^ * total item failure rate 


total item operating time per year, determined by 
the ACF. 


Tables IV and V tabulate the calculated results for X^ = NX and 
Table V contains the results for F. 


4.3.3 DOWNTIME OR RECOVERY TIME . Downtime is the time it takes to restore 
system operation after an in-service failure or malfunction. The 
cumulative downtime per year for a given item is determined by 


F x t 


where: F = number of item failures per year 


t f * downtime or recovery time per failure. 


The calculated results of the cumulative downtime for the subsystems 
and items are summarized in Table V. 


> 
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TABU V-l. PREDICTED ITEM EATU1R U ’ COUNT AND DOWNTIME (ACE - O.lo) 



PLANT SIZE = 0. 

5MWe 

PLANT STZE = 1. 

OMWe 

PLANT SIZE - 10. 

OMWe 

SUBSYSTEM/ ITEM 

■H' 

t 

X. 

F 


X 

F 

nwm 

X 

n 

Fx t 


iHfnHi 

r 

I 



1 



1 


- L 

1 . COLLECTOR 








u 


mm 


A. Structure 

2838 

24.0 

3.3 



6.6 

0.02 

■ 

66.6 


4.8 

B. Drive Controls 

2838 

2.0 

548.1 

1.5b 


109b. 2 

3.11 


10962.0 


b2.2 

C. Mirrored Panel Set 

2838 

8.0 

112.5 

Blfl 


225.0 

0.64 


2250.0 


51.2 

D. Receiver/Vapor Pipe 

2838 

5.5 

150.3 

mJLvm 

2.37 

300.6 



3006.0 

8.5 

46.8 

COI .LECTOR TOTAL 




2.32 

8.29 


4.62 

lb. 50 


4o.2 

165.0 

2 . POWER CONVERSION 












A. Engine and Access. 

2838 

4.0 

697.5 

1.98 

7.92 

1395.0 

3.9C 

15.84 

13950.0 

39.6 

158.4 

B. Alternator & Acces. 

2838 

4.0 

165.2 


1.88 

330.3 

0.94 

3.76 

3303.0 

9.4 

37.6 

POWER CONV. TOTAL 




2.45 

9.80 


4.9C 

19.60 


49.0 

196.0 

3. 'J'RANSPORT/DISTRIB. 












A. Main Switchboard 

2838 

0.5 

86.4 

0.25 

0.13 

164.2 


0.24 

1563.8 

4.4 

2.2 

B. HV Transformer 

2838 

12.0 

5.7 

0.02 

0.24 

5.7 


0.24 

57.8 

0.2 

2.4 

C. HV Switch 

2838 

1.5 

1.6 

0.00 

0.00 

1.6 


0.00 

16.0 

0.0 

0.0 

TRANS /DTSTR. TOTAL 




0.27 



0.4'. 

0.48 


4.6 

A . b 

4 . COMPUTER CONTROL 












A. Computer & Access. 

2838 

1.0 

155.5 

0.44 

0.44 

155.5 

0.44 

0.44 

155.5 

0.4 

0.4 

. B. Fan/Air Cond. 

2838 

1.0 

216.0 

0.61 

0.61 

216.0 

0.61 

0.61 

216.0 

0.6 

0.6 

COMPUTER CONTROL TOTAL 




1.05 

1.05 


1.0' 

1.05 


1.0 

1.0 

SYSTEM TOTAL 




6.09 

19.51 


11 .Of 

37.63 


100.8 

'66.6 


Legend: 

= item operating hour in hours per year. 


t = estimated item recovery time per failure 
= downtime in hours. 

X. = item base failure rate in .failures per million hours 
for a given plant size and ACF (See Table IV) = NX 

F = (X ) ( t Q ) = failure quantity over the duration. 

l’xt = item downtime over the t duration. 
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1 . COLLECTOR 

A. Structure 

B. Drive Controls 

C. Mirrored Panel Set 

D. Keeeiver/Vapor Pipe 

COLLECTOR TOTAL 

2 . POWER CONVERSION 

A. Engine & Access. 

U. Alternator & Access. 

POWER CONV. TOTAL 

3. TRANSPORT/DISTRIB. 

A. Main Switchboard 

B. 11V Transformer 

C. HV Switch 



TRANS/ DI SIR. TOTAL 

4. COMPUTER CONTROL 

A. Computer & Access. 

B. Fan/ Air Cond. 


5 . STORAGE 

A. Fan/Air Cond. 

B. AC to DC Conv. 

C. Battery Bank 

D. DC to AC Inverter 
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TABLE V-3 . PREDICTED ITEM FAILURE COUNT AND DOWNTIME (ACI- - 0.70) 


PLANT SIZE = 0 . 5MWe PLANT SIZE = l.OMWc 


PLANT SIZE - lO.UMWe 


SUBSYSTEM/ ITEM 


1 . COLLECTOR 

A. Structure 

B. Drive Controls 

C. Mirrored Panel Set 

D. Receiver/Vapor Pipe 

COLLECTOR TOTAL 

2 . POWER CONVERSION 

A. Engine and Access. 

B. Alternator & Acces. 

POWER CONV. TOTAL 

3. TRANSPORT/ DISTRIB. 

A. Main Switchboard 

B. 11V Transformer 

C. 11V Switch 

TRANS/ DISTR. TOTAL 

4. COMPUTER CONTROL 

A. Computer & Access. 

B. Fan/Air Cond. 

COMPUTER CONT. TOTAL 

5. STORAGE 

A. Fan/ Air Cond. 

B. AC to DC Conv. 

C. Battery Bank 

D. DC to AC Inverter 

STORACE TOTAL 


SYSTEM TOTAL 


7.8 

1278.9 

262.5 

350.7 


1627.5 

385.4 


198.7 

5.7 


155.5 

216.0 


216.0 

101.0 

278.8 

100.0 




0.96 151.7 
14.18 24969.0 
11.00 5125.0 
10.67 6847.0 


36.08 31775.0 

8.56 7523.5 


3637.4 

114.0 

32.0 


216.0 


4320.0 

2020.0 

5576.0 

2000.0 


105.2 374.1 


90.2 360.8 

21.4 85.6 


111.6 446.4 





896.4 


























4.3.4 CUMULATIVE-MEAN- TIME- 3E TWEEN-FAILURE . The CMTBF is defined by 

the expression 

CMTBF =* T o 


where: 


total time a given system delivers 
rated power output per year as 
determined y the ACF, and 


IF * cumulative number of failures that 
are expected to occur during system 
operation for a year. 

The CMTBF for the nine systems are tabulated in Table I. 


4.3.5 CUMULATIVE MEAN- TIME - 
the expression 


£Y. The CMTTR is defined by 


CMTTR 


Z£Fxt r ) • 
IF 


The CMTTR for the nine systems are tabulated in Table I. 

4.3.6 SYSTEM AVAILABILITY. The system availability or up time 
probability is defined by the expression 

A - CMTBF . 

CMTBF + CMTTR 

The system availability results are tabulated in Table I. 


f 
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APPENDIX J 


r 



STUDY OF RECEIVER INSULATION MATERIALS 


Table J-l presents some of the physical and cost data for various thermal 
insulations investigated for SPS. The table was limited to high temperature, 
low thermal conductivity, low cost and lightweight insulations. These in- 
sulators are typically ceramic blankets made from high purity alumina- silica 
fibers, and the data were obtained from several vendors. The insulation cost 
is listed in terms of dollars per board foot (1" x 12" x 12"). This cost is 
the "best" quotation price obtained for a purchase of 5,000 (or more) board 
feet of insulation. This cost can be reduced by purchasing insulation in 
greater quantities and by purchasing directly from the manufacturers. 

Thermal conductivities of the insulators are presented as a function of mean 
insulation temperature in Figure J-l. The thermal conductivities for the 
ceramic blankets at the mean temperature (~375°C) of the receiver subsystem 
are approximately 0.06 to 0.08 W/m-°C. The thermal conductivity of the Johns- 
Manville Min-K is approximately one-half that of the other insulators, jw 
ever, it is more appropriate to compare thermal insulators on a "performance- 
cost" basis. Insulation performance-cost is simply the product of the 
thermal conductivity and the cost per volume. Figure J-2 presents this factor 
for several candidate materials as a function of the mean insulation tempera- 
ture. The "Cerawool"* and "Durablanket" insulations are the least expensive 
ceramic blankets in the 8#/ft3 (128 kg/ttP) densities. The Cerawool blanket 
has been selected for the baseline since it has slightly lower conductivity. 
Although less expensive, the lower density insulations such as the "Duraback" 
blanket are not recommended since they are easily compressed to a less effec- 
tive thickness. The performance-cost of the Cerawool blanket is approximate- 
ly 35$-W/nP-°C. Although it has a very low thermal conductivity, the per- 
formance-cost of the Min-K is about 1350$W/m" + -°C, which makes it prohibitive- 
ly expensive. 





i 


* Manufacturers of these trade-name insulations are listed in Table J-l. 

I 
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TABLE J-l: INSULATION DATA 

MAXIMUM 

CliQfICAL SERVICE 



INSULATION 

COMPOSITION 

TEMPERATURE 


DENSITY 

COST 

1 

Cerawool A 

40% A1 2 0 3 , 

51% 

S i°2 

1600° F 

(871°C) 

8# /ft 

3 (128 kg/in 3 

1.08 $/ BDf t 

2 

A 

Cerablanket 

47% A1 2 0 3 , 

53% 

S i°2 

2400 

(1316) 

8 

(128) 

1.42 

3 

Cerafelt A 

46% A1 2 0 3 , 

54% 

S i°2 

2400 

(1316) 

8 

(128) 

1.70 

4 

Cerachrome A 

41% A1 2 0 3 , 

55% 

S i°2 

2600 

(1427) 

8 

(128) 

3.04 

5 

D 

Durablanket 

48% A1 2 0 3 , 

52% 

S i°2 

2300 

(1260) 

8 

(128) 

1.14 

6 

Duraback B 

40% A1 2 0 3 , 

57% 

S i°2 

1700 

( 927) 

4 

( 64) 

0.43 

! 7 

g 

Lo- Con Blanket 

52% A1 2 0 3 , 

48% 

S i°2 

2300 

(1260) 

6 

( 96) 

1.08 

8 

Mineral Wool C 




1000 

( 538) 

8 

(128) 

0. 15 

9 

Krawool^ 

45% A1 2 0 3 , 

52% 

S i°2 

2300 

(1260) 

8 

(128) 

1.65 

10 

Fiberglass TIW-I E 




1000 

( 538) 

- 

— 

0.09 

n 

Min-K* 




1800 

( 982) 

16 

(256) 

84.00 


A. Johns-Manville 

B. Carborundum Company 

C. Holmes Ltd. 

D. Babcock & Wilcox Company 

E. Owens-Corning Fiberglass Corporation 



| CONDUCTIVITY, k, w/m- 



NUMBERS REFER TO 
MATERIALS IDENTIFIED 
IN TABLE J-I 


2500 ; 


FIGURE J-l. THERMAL CONDUCTIVITY OF 
INSULATION MATERIALS 


IS 

QUALITY 










NUMBERS REFER TO 
MATERIALS IDENTIFIED 
IN TABLE J-I 


1000 (OF) 1500 
MEAN INSULATION TEMPERATURE 

00 600 (°C) 800 100( 


FIGURE J-2. PERFORMANCE-COST OF 
INSULATION MATERIALS 
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BUFFER STORAGE MATERIALS AND COMPARISON 












APPENDIX K 

BUFFER STORAGE MATERIALS AND COMPARISON 


1. MATERIALS 

3uffer storage is composed of a material capable of absorbing and releasing 
large amounts of heat over the proper temperature range. Storage capacity can 
be in the form of sensible heat, latent heat, heat of reaction, or various 
combinations thereof. The work of Schroder at the Philips Aachen Lab 
(Germany) demonstrates that the most practical solution at the present time 
for Stirling engine applications is a combination of sensible energy of the 
solid and liquid phases plus latent heat of fusion of certain high-melting- 
point salts (Reference K-l). This is verified to some extent by a Sandia 
study (Reference K-2) in which some additional compounds were identified 
which may have some cost advantages. 






Investigation of the basic laws of chemistry shows that the highest sensible 
heat capacities and heats of fusion per unit weight and volume will ob- 
tained in substances having a close molecular packing of small and light 
elements. (Reference K-3) . In addition to high thermal capacity, a storage 
material should satisfy several other criteria : 


(1) The melting point must be within the operating temperature 
range of the engine, which should be as high as possible for 
maximum power plant efficiency (i.e., ~800°C for the exist- 
ing USS Stirling engines) . 




(2) The material must be chemically stable and non-corrosive or 
have additives to eliminate corrosion. 

(3) Low cost, especially for large-scale applications. 

(4) The vapor pressure should be low at the maximum temperatures. 

(5) The material should not be highly flammable or toxic. 



The first criterion sets a fusion point of the buffer store in the range of 
roughly 750-S50°C (1380-1560°F) . It is usually desirable to run a Stirling 
engine at a high constant heater head temperature as much as possible since 
the engine operates less efficiently at lower heat temperatures. Since this 
can occur only at the stores fusion point if the engine is running exclu- 
sively from the buffer heat storage, it is desirable to have a large heat of 
fusion. The lower temperature limit at which the buffer storage is effective 
is dictated by the heat transfer characteristics of the medium which trans- 
ports the energy from the store to the engine. 

Studies documented in References K-l and K-5 have shown that the fluorides 
of lithium, sodium, and magnesium have higher total heat capacities than any 


3 




K-l 


other practical materials. However, the high melting point of most pure 
fluorides (except LiF and BeF2) makes it necessary to employ lower melting 
point eutectic mixtures to obtain the desired performance over the selected 
operating temperature range. Other properties of fluoride-type materials 
are also as good as, or better than, most other potential heat store candi- 
dates, particularly for the lithium- based compounds. Unfortunately, the 
current price of lithium compounds appears to make them non-competitive on 
a cost basis for large-scale application, and beryllium compounds are ex- 
cluded on the basis of cost and toxicity. Other potential candidates, 
cirticularly from a cost standpoint, are basically in the chloride group 
(Reference X-2) . The following high temperature heat store salts were 
selected for comparison from a much larger list of candidates. 

Fluorides : LiF, and eutectics of LiF/MgF2, NaF/MgF 2 , or LiF/NaF/MgF 2 

Chlorides : NaCl, NaCl/CaCl 2 eutectic, NaCO^/KCl eutectic and CaCl^ 

Carbonates : K 2 C 03 /Na 2 C 03 eutectic 

Fluorides . Table K-l lists some of the properties of candidate fluorides. 

LiF has the highest heat of fusion of all the materials on a weight basis but 
contains a large amount of scarce and expensive lithium. An eutectic mix- 
ture of 67" LiF/33% MgF 2 (mole percent) contains only 12.1 percent lithium, 
costs about one-fourth that of LiF, and has a volumetric heat of fusion 
even higher than LiF (or even LiH) . The third mixture shown in the table 
contains no lithium and is one of the more cost effective candidates, as 
well as having other favorable properties as outlined below. The low pro- 
duction cost of the NaF/MgF2 material is due to the fact that fertilizer 
manufacturers have been throwing away fluoride-bearing compounds (mostly 
in the form of ^SiFg) for a number of years as waste products. Philips 
has obtained cost quotes for the production of material No. 3 of $0. 13/kg 
in large quantities from LaPorte Company in England. This quote includes 
the recovery of fluorine and combining it with common sodium and magnesium 
(or other elements as desired) . The production cost of this material in the 
U. S. is also projected to be about $0. 15/kg.*, provided fluorine can be 
obtained from ^SiFg, which is now an undesirable waste product and is 
usually dumped in a land fill. Today the price of NaF or NaF/MgF 2 is about 
$0. 70/kg* because the source of the flourine is expensive hydrofluoric acid, 
and production is relatively small. The fourth material is typical of a 
relatively low melt temperature fluoride eutectic, and is included only for 
reference since the cost will likely be too high and the melting point too 
low compared to material No. 3. 


*Personal communication with Dr. J. Eichelberger of the Pennwalt Corp. 
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TABLE K-l. REPRESENTATIVE PROPERTIES OF CANDIDATE HEAT STORE MATERIALS 
(Sources: References K-l, K-2, K-5, K-7) 


i 


Composition, 
No. Mole Z 

1 LiF 

2 67 LiF/33 MgF 2 

3 75 NaF/25 MgF 2 

4 46LiF/44 NaF/ 
10 MgF 2 


Melting Allf, Heat Heat Capacity; Density at Melt, Approx. 

Point oi Fusion 450-850°C __ g/ cm Cost 

°C (°F) kW-sec/kg(BTU/lb) kW- sec /kg (BTU/lb) Solid 1 Liquid" $/ kg 


B48 (1558) 930 

746 (1375) 910 

832 (1530) 636 


930 (l * 

(400) 1885 

910 

(390) 1640 

636 

(273) 1275 

686 

(295) 1260 


1.83 
2.63 2.31 
2.69 2.19 
2.61 2.10 


20-22 
/. on ) 


Comments 

Too expensive 

Too expensive 

Final candid. 

Probably too 
expensive, 
melts too low 


5 NaCl 


80 ^ ) 


(1472) 479 


6 48 NaCl/ 52 CaCl 2 500 (932) 281 

1 50 Na 2 C0 3 /50KCl 588 (1090) 268 

8 CaCl 2 772 (1422) 255 


(206) 1020 (439) 

(121) 720 (309) 
(115) 830 (357) 
(110) 650 (280) 


9 44 K2CO3/56 Not evaluated, 

710 (1310) 159 ( 68) judged to be 

Na 2 CQ 3 low 

(1) Clven as 1046 kW-sec/kg (450 BTU/lb) in some references. 

(2) Reference K-6; unknown quantity and impurity level. 

(3) Projected production costs (see text). 

(4) Reference K-7; other references give melt point as high as 808 C. 

(5) Reference K-3; unknown quantity and impurity level. See text for comments. 


Na 2 C0 3 


710 (1310) 159 


1.96 1.55 See Text 

G.073* 5 * 

o.n< 5 > 

2.09? 2.08 See Text 


Possible candid, 
but low P 

Melts too low 

Melts too low 

Possible can- 
didate but 
low heat 
capacity 


Heat capacity 
too low 



Additional advantages of the f luoride-type compounds are as follows. First, 
the densities of the materials near the melt are considerably higher than 
chlorides or other competing compounds. This is a significant advantage for 
packaging. Also, the vapor pressure is low (less than 1 Torr up to 900°C) , 
they have high chemical stability, and a thermal conductivity which is higher 
than most other types of storage materials (solid conductivities range from 
about 4 to 8 W/m-°C). However, the conductivity is still relatively low, 
which necessitates a package with a large surface area and a low heat trans- 
fer rate in order to avoid excessive temperature gradients in the salt 
during the periods when it is solidified. Measurements have shown that the 
expansion coefficients are nearly linear and there is no phase change in 
the solid state. The densities of the individual solid fluorides differ 
from one another and from the densities of the eutectic melts by less than 
20 percent. No supercooling is evident, so that there is little tendency 
for the phases to separate by sedimentation. Because of their high chemical 
stabilities, the fluorides are rather inert and show little corrosive action 
on stainless steel. When a small amount of aluminum is added to the fluoride 
melt as a corrosion inhibitor, no corrosion has been observed on 18-8 type 
stainless steel after well over 1000 hours of operation at 850°C in the 
Philips Eindhoven (Holland) Laboratory. An inert atmosphere is required at 
these conditions. 


Chlorides . Materials 5 through 8 in Table K-l list some typical chlorides. 
The two eutectic compositions (Nos. 6 and 7) have far too low a melting 
(fusion) temperature, and are not considered further. Both sodium chloride 
(No. 5) and calcium chloride (No. 8) have melting points in the proper range. 
They have heat capacities somewhat below material No. 3, especially for 
calcium chloride. Other comparative disadvantages include: lower densities, 

apparently lower thermal conductivity (although good da^a are lacking), and 
the lack of practical experience in their use as a high temperature buffer 
store. 


There are potential corrosion problems with chlorides unless pure, anhydrous 
grades are used in an inert atmosphere. An initial examination indicates 
that USP grade NaCl should be used pending corrosion studies of less pure 
grades. No water should be present, although if some occurs, it may not be 
too great a problem to drive it off by heating. Anhydrous grade CaC^ must 
be used, again, pending further tests since wet CaCl 2 is expected to have 
corrosion problems with stainless steel containers. Heating is not a good 
way to dry CaC^ since oxides tend to form which are also quite corrosive. 


The cost of NaCl and CaC^, derived from Reference K-2 is between 7 and 9q/kg. 
This is as low a cost as can be expected for any chemical compound, thus these 
materials may have a cost advantage for use as heat stores. They will main- 
tain this advantage over, say, material No. 3 as long as they cost less than 
approximately % to 1/3 that of NaF and MgF 2 . However, it appears that 7-9q/ 
kg is too optimistic for the high purity grades needed to assure that cor- 
rosion is not a problem. The costs in 1976 for the Los Angeles area are 
as follows: 
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Material 

Grade 

Total Amounts 

Containers 

Cost 

L 

NaCl 

l T SP 

Truckload 
(24,000 lbs.) 

350 lb. drums 

30.37/kg 

1 

CaCl 0 

Anhydrous 

Truckload 
(24,000 lbs.) 

350 lb. drums 

$0.64/ kg 


Both materials are available in less pure grades (technical grades) in 
freightcar lots or greater at much lower prices, but these grades are not 
recommended until the corrosive nature is actually tested at the proper 
temperature range and/or pre-treatment techniques are investigated. 

Carbonates . The ninth material identified is a eutectic of potassium car- 
bonate and sodium carbonate. A very low heat of fusion makes this material 
undesirable for buffer store applications compared to the other leading 
candidates • 

Material Comparisons . Figure K-l shows a comparison plot of the heat store 
capacity of several typical materials over the selected temperature range. 

In most cases the sensible energy and the latent heat of fusion are roughly 
equal, and the heat store would have to be correspondingly larger if sensible 
energy were not used. As discussed previously, the recommended temperature 
range is dictated by the heat transport capability of the heat transfer 
fluid on the low side, and the temperature limits of the Stirling engine on 
the upper end. 


Figure K-2 presents the cost of heat store capacity for selected materials, 
neglecting container and all other non-heat-st.ore-material costs. The results 
for mixtures No. 6 and 7 are based on the costs derived from Reference K-2 
which appear to be quite optimistic. However, note that materials No. 6 and 
7 are not candidate materials because of their low fusion point, but the plots 
demonstrate the potential cost/storage energy relationship of chlorides if 
they could be purchased in the 7-llc/kg range. Material No. 3 is quite com- 
petitive on a projected production cost basis of 15c/kg., especially consi- 
dering its other advantages. One of the best materials from strictly a total 
heat capacity standpoint, LiF/MgF-,, is very poor on a cost comparison be- 
cause of the $4.98 /leg figure. 

In some cases it may be advantageous to change the fusion point of the 
storage material to a slightly lower value than exhibited by a single, pure 
compound. This is achievable with a non-eutectic mixture, and range of 
fusion temperatures is exhibited, rather than a single value. There are 
some well-known problems associated with the precipitation of different 
compounds at different temperatures for some mixtures, but their suitability 
should be investigated further. 


For example, the NaCl/CaCl 0 phase diagram shows that a mixture of 80 NaCl/20 
CaCl 2 (mass fraction) will begin to precipitate almost pure NaCl (actually, 
a solid solution of CaC^ in NaCl) at ~760°C (1400°F). As this separates 
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FIGURE K-2. COST OF HEAT STORE CAPACITY FOR 
SELECTED MATERIALS 


from the remaining liquid, the liquid will become richer in CaC^* The solid 
solution will concinue Co separate until the peritectic temperature and com- 
position are reached (~600 J C and mole percent NaCl) . At this tempera- 
ture che liquid will react with Che precipitated solid solution to form the 
compound iNaCl-CaC^ until che melt is completely solidified. Thus, fusion 
starts at ~760°C and is completed at ~600°C for this particular mixture. 

In summary, the leading candidates for the high temperature energy storage 
material are in current order of preference: 

1. 75 NaF/25 MgF-, 

2. NaCl 
3 • CaC 1 2 

4- Non-eutectic mixtures of NaCl and CaCl 2 

Considerably more systems and cost analyses will be required to verify these 
initial choices 

2. ENCAPSULATION 

The buffer storage system must efficiently transfer heat from the heat store 
material to the heat transfer or working fluid (and from the working fluid 
into the mixture during charging). N. V. Philips has investigated this 
problem in terms of various encapsulation techniques and determined chat a 
good design is one that has a large number of long, slender, thin-wall stain- 
less steel containers which are filled with the score material and sealed. 
These containers are in turn enclosed by liquid sodium within the receiver 
pool. A photograph of a partial assembly of a laboratory heat store showing 
Che individual containers, or canisters, is presented in Figure K-3. When 
the mixture is liquid it occupies almost the entire volume in the canister 
except for a small volume of inert argc r ’ gas- When the mixture solidifies 
there is a slight argon over-pressure so that the container will not collapse. 
Figure K-4 shows a cross-section of a canister filled with solidified 
lithium fluoride, and Figure K-5 shows the external view of the tank which 
holds the canisters. The tank insulation is partially removed to show the 
external configuration and the many temperature and pressure instrumentation 
leads . 

During charging, energy from the solar heated sodium pool melts the 
eutectic salt. Whenever solar power is not received (and the salt is 
charged), the latent heat of fusion is released into the pool thereby vapor- 
izing sodium to continue engine operation. 

3. COST OF BUFFER STORAGE 

The cost of buffer storage is obviously a function of the material which is 
selected -- not only because of the store material cost, but also because of 
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fr- 
its heat capacity and density characteristics. For example, Table K-l shows 
that the heat capacity of NaF/MgF 2 is 25 percent better than NaCl, and the 
density of the solid at the melt point is 37 percent higher (41 percent 
higher for the liquid phase). Therefore, NaCl will require a 70 percent 
increase in store volume compared to the baseline material for the same 
energy storage (assuming other conditions are equal). It is obvious that 
the baseline material can be a factor of 2 or 3 more expensive on a weight 
basis than common salt and still be competitive. Another important cost 
factor is the design and fabrication of the canisters. Low density, lower 
heat capacity materials require larger or additional containers, thus in- 
creased containerisation costs. 
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FIGURE K-3 


PARTIAL ASSEMBLY OF PHILIPS HEAT 
STORE SHOWING CONTAINERS FILLED 
WITH STORE MATERIALS 


FIGURE K-4 


CUT-AWAY OF A CANISTER SHOWING 
SOLIDIFIED HEAT STORE MATERIAL 
(LITHIUM FLUORIDE) 










FIGURE K-5. PHOTOGRAPH OF PHILIPS HEAT STORE TANK 
AND PARTIAL INSULATION COVERING. 
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1. 0 INTRODUCTION 

This Technical Report presents the development of an improved transient 
thermal model for the baseline sodium receiver and thermal transport sub- 
system. An in-depth transient analysis of the receiver-thermal transport 
subsystem is necessary in view of the inherent non-steady nature of the 
solar conversion process. The objective is to characterize the behavior 
of the subsystem for all operational modes encountered in a solar application. 
Thus, the preliminary thermal model presented in Reference 1 and an unpub- 
lished transport model have been reformulated into a new detailed trar jxent 
thermal model. This improved transient analysis incorporates the physical 
and thermal characteristics of the latest baseline design. Subsystem re- 
sponses to normal and inclement operating conditions are presented for the 
system using the USS P-75 Stirling engine operating at 1800 rpm. 

2.0 METHODS OF ANALYSIS 

A schematic of the thermal model is presented in Figure 1. As shown, the 
energy balance considers the incoming solar power, the thermal power to the 
engine, and the thermal losses from the subsystem. For computational con- 
venience, conduction losses from the vapor pipe were distributed equally 
between receiver and engine heater head. Pressure drops throughout the 
entire subsystem have also been included. 




2. 1 MODEL IMPROVEMENTS 

The following are the major improvements made to the transient thermal model. 

- Explicit integration with respect to the time variable. 

- Separate differential equations for the pool and head temperatures. 

- Deletion of the sodium vapor thermal capacitance since it is 
negligible. 

- Inclusion of vapor transport pressure losses. 

- Consideration of various engine, receiver lid, and vapor valve 
control options. 

- Refinement of the solar input heating function to include a clear 



- Refinement of the engine output function. 

- Elimination of the minor numerical instabilities encountered during 
the development of the model. 

- Streamlined computational sequence. 

- Refined input/output formats. 

2.2 MODEL DEVELOPMENT 

The nomenclature and units for the following model are presented in Appendix 
I. It is assumed that the sodium vapor is saturated (vapor and liquid are 
in quasi-equilibrium), the vapor volume is approximately constant and that 
receiver orientation effects are negligible. 

Basic energy balance: 

Solar 3 Q RR + SoNV + SoOL + SoNDH + q head + q eng 


Basic differential equations: 

The following two differential equations are solved in the computer code 
using a library subroutine for numerical solution of simultaneous first- 
order ordinary differential equations with automatic step change. 


d T. 


POOL 

dt 


(Q 


SOLAR 


Sr ‘ SoNV 


SoNDP " ^ TRANS ^ 


c sodium + Structure + c salt 


d T. 


head 

dt 


(Q TRANS ' SoNDH ' Q ENG^ 
C HEAD 


Solar input heating function: 

The following equation generates the solar input power to the receiver. 

Solar 3 £ drop Sjish ^dsh I o exp( " K ABs f abs * 


>3 



where: 


ABS 


I 

cosL cos; cos [tt(1 * t/720)' + sinL sin^ 


This solar heating function is based upon the 1976 15 minute Barstow 
insolation data. A typical uncloud> day was curve-fit to determine 
the clear sky absorption term. With the present code, the solar input 
heating can be characterized with any combination of clear, hazy and 
inclement conditions as a function of solar time. 

Heat losses: 

Heat loss from the subsystem is characterized by the following equations: 
Q RR " (ffei V (T W ’ T SKY ^ f LID 





Q CONV 

q cont>p 

q condh 


< h V (T V-W f LID 

~ G POOL + G LID (1 " f LID^ - (T POOL 
- G HEAD + G ENG (1 " f ENG^ (T HEAD 


T 

AMB 

T 

AMB 


) 

) 


where: 



W 


t pool + (q solar ’ q rr q con\P /g wall 

1 /2 V 

+ - (Q SOLAR “ G RR ‘ Q CONV^ ^Sla- 


The effect of using sky, ambient or zero sink temperature for the 
re-radiation heat loss term is negligible for the baseline subsystem 
(only a ± 25W difference). 
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Power transported: 

The thermal power transferred from the receiver pool to the engine 
heater head is given by the following equation. 



and 


( Vc 


HEAD 


'POOL 


P POOL 'SlJCT g 
* 
c 


(ZRgTpo OL /M) 
“ P REF “P 
" P REF «P 


1/2 j 

“*7 

U 


r 




T 

HEAD 

T 

REF 

t pool 

T R£F 
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The sodium heat of vaporization is approximated by the following equation: 



&H - 4.9815 x 10 6 - 1022. 7T 

V 



W - sec 

ES 


Figure 2 compares equation (1) with the data given in Reference 2. Figures 
3 and 4 present the sodium vapor pressure (P^) and ratio of specific heats 
(y) as functions of temperature, respectively. The following equation is 
used to approximate the sodium gas constant. 


- 296.3 + 65.5 cos [tt/2 ( - 0.5)] 


ZR W - sec 
M ~ k * " K 


( 1 ) 


( 2 ) 


This equation is compared with the exact value in Figure 5. 


| 



Engine input power: 

The engine power parameters are based upon the analyses of Reference 4 
and are presented below. 


'ENG 


£ 100 


100 / (pr) head 


(PR) AMB f ENG 


where: 

hea d " 


- 273) 0,2265 - 3.545 

(Thead - 273 >°' 2U9 - 
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2.0224 - 0.00355 T 


(PR) 


COOL 


4MB 1.7978 - 0.00277 T 


COOL 


and where: 

T COOL “ T AMB + ^ENG “ Q OUT )/G COOL 


and where: r 
Q 0UT 


i T 


Q o + 2 Q 100 V 100 


* hi 


NG 


* [ "head ' 273 >°‘ 2265 ' 3 


.545 j 


[*■ 


0224 - 0.00355 T 


COOL 


n Q 

ENG y ENG 


The Receiver-Thermal Transport Model presented above has been programmed for 
operation on the Aeronutronic Time-Sharing System (TSS) . 


3.0 RESULTS 

Detailed transient results are presented in this section for the baseline 
receiver-thermal transport subsystem with the Stirling engine operating at 
the baseline 1800 rpm condition. Table I presents some of the physical and 
thermal data of the subsystem used for the present analyses. Both normal 
and inclement operating conditions have been analyzed. 

The following terms are used in this section and are defined as follows: 
"Idle", "rated" and "equilibrium" power refer to the Stirling engine shaft 
output power level. 

"Idle " power is the low power necessary to stay "on-line". 
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"Rated " power refers to the power required per engine for 
the system to generate rated electrical power (1 MW e ) at 
rated solar conditions (I, ■ 800W/m“, T tvrD ■ 44. 6°C). 

" Equilibrium " power is the power level required by the 
engine controller to maintain the steady-state head 
temperature. 

" Buffer storage " refers to thermal storage added to the 
subsystem, and is therefore in addition to the inherent 
thermal capacitance of the subsystem, 

• NORMAL START-UP 

Figure 6 presents the subsystem temperature response and engine output power 
during a normal start-up condition. The solar input power profile is based 
on the 15 minute Barstow insolation data for a typical uncloudy day and in- 
cludes the effects of both concentrator size and efficiency, thus representing 
the total power entering the receiver. This solar profile has been used for 
all start-up cases to facilitate direct comparison between them. In the 
morning, after normal operation from the previous day and nighttime cool down, 
the receiver sodium pool temperature is predicted to be approximately 525°C 
and the engine head at approximately ambient temperature. The receiver lid 
and vapor valve are opened prior to focusing on the sun. Once the sun is 
"on", the pool and head temperatures rapidly increase towards the steady-state 
operating level. After about 3 minutes the engine head temperature has reached 
the predetermined start temperature (450°C for the present analysis). The 
engine is then started and operated at a low or "idle" power level consistent 
with a previously selected helium mean working pressure (F) for idle conditions 
(2.5 MPa for the present analysis). When the engine is started, an Inflection 
of the receiver pool and engine head temperature slopes is observed. This 
occurs because some of the heat, which would normally Increase the temperatures 
is now being removed by the engine. Since the net power to the receiver is 
increasing, the two temperatures continue to Increase while the engine operates 
at "idle" power. After approximately 8 minutes the engine head reaches the 
steady-state operating temperature (800°C) and the engine head temperature 
control mode is activated. At this point the temperature controller moves 
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Che power valve to bring the engine P up (and hence the torque level) to an 
"equilibrium' 1 power level. The power valve is then modulated to maintain 
the engine head temperature within the control band as the solar flux varies 
during the day. Reference 5 presents a detailed discussion of the solar- 
Stirling engine control scheme. 

• NORMAL SHUTDOWN 

Figure 7 presents the subsys :em temperature response during normal shutdown 
in the evening. For this case, the engine is assumed to be operating at the 
steady-state temperature and at any "equilibrium" power level above "idle" 
until the incoming solar power drops to zero. At this time the receiver lid 
is closed and the engine continues Co produce power from the inherent thermal 
capacitance of the receiver. As energy is removed, the engine head tempera- 
ture drops below the steady-state level. The engine is then operated at the 
previously mentioned "idle" power level. After 15 minutes the head tempera- 
ture reaches the lower operational limit and the engine and valve are shut 
off. The subsystem temperatures will then slowly soak out to the ambient 
temperature. Another mode for normal evening shutdown is presented later. 

• NORMAL START-UP WITH BUFFER STORAGE 

Start-up from the previous day with the subsystem containing 100 kg of buffer 
storage is presented in Figure 8. Three different eutectic salt melting temp- 
erature regions have been investigated in the present analysis (i. e. , 790-800°C, 
802°C, 803. 5°C). Approximately four minutes after the sun is focused upon 
the receiver the engine head temperature reaches the lower operational limit 
and the engine is started. An additional twelve minutes is required before 
the engine is controlled by head temperature for the 802°C and 803. 5°C salt 
melt temperatures. The time required to melt the 100 kg of buffer storage 
for the 802°C and 803. 5°C salt melt temperature cases is 26 and 62 minutes, 
respectively. For cases where the salt melt temperature is between 790°C 
and 800°C, approximately eleven minutes of "idle" engine operation is required 
while the 100 kg of buffer storage is melting. The engine head and receiver 
pool temperatures rise once again after the eutectic salt is entirely melted. 
When the head temperature reaches the steady-state operating temperatures of 
800°C the "idle" control is shifted to head temperature control. Comparing 
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Figures 6 and 8 shows that the temperatures rise more slowly with the buffer 
storage, as is expected. (Note scale change.) For example, three minutes 
are required for engine 3 tart-up for the baseline subsystem and four minutes 
are required for the subsystem with 100 kg of buffer storage. 


• EXTENDED "BAD" WEATHER 


Figure 9 presents the subsystem temperature response during normal nighttime 
and extended "bad" weather conditions. For this case the engine and valve 
are assumed to be shut off as soon as the solar power becomes zero. With 
the receiver lid open the pool temperature drops to the sodium freeze temp- 
erature (98°C) in about sixteen hours. However, if the receiver lid remains 
closed, the pool temperature is predicted to stay above freezing in excess 
of three days. It is interesting to note, based upon the site data, that 
the sodium pool would not have frozen at any time during 1976 at Barstow. 

In the absence of additional engine coolant pumping after shutdown, the head 
is predicted to reach ambient temperature in about a day. The head temperature 
is independent of the pool temperature response since the closed vapor valve 
effectively uncouples the head and pool. As shown in Figure 7, the pool and 
head temperatures at the time of engine shut-off are 545°C and 450°C, respec- 
tively. By using these two temperatures as initial conditions, the soak-out 
(or non-operating cool down) temperature response for the normal nighttime 
shutdown presented in Figure 7 can easily be obtained from Figure 9. 

• START-UP WITH FROZEN POOL 

Start-up with a frozen sodium pool can occur after extended "bad" weather 
conditions or during the initial start-up. The subsystem temperature start- 
up response for this condition is presented in Figure 10. Engine start-up 
and switch to head temperature control occur after about 10 and 15 minutes, 
respectively. As shown, less than one minute is required to melt the sodium 
pool. As the sodium melts, a liquid layer forms between the heated receiver 
wall and the still-frozen sodium. This liquid layer rises via capillary 
action exposing solid sodium to the heated surface. No problems are ex- 
pected during the melt period based upon the extensive work by the N. V. 

Philips Company (Reference 6). However, if any problems do occur when the 
sodium is melted directly by the sun, nichrome heater wires can be placed 

>10 

\ 



between the outer receiver wall and the insulation, and current passed through 
the wires to slowly melt the entire sodium pool before the receiver is exposed 
to the sun. This method requires approximately 8. 25 kW-hr of electrical 
energy to raise the receiver temperature from ambient to a level above the 
sodium melt temperature. 

• CLOUD PASSAGE 

Figures 11 through 14 present the subsystem temperature response and engine 
output power during cloud passage conditions. The same solar input power 
profile (based upon the Barstow insolation data) has been used for all cloud 
passage cases to facilitate direct comparison between the various cases. 
Steady-state temperatures and "equilibrium” power conditions are assumed to 
exist prior to the onset of cloud passage. Cloud passage was assumed to 
start at 11:00 a.m. for all cases presented. 

• 15 MINUTE CLOUD PASSAGE - BASELINE SUBSYSTEM 

Figure 11 presents the baseline subsystem temperature response to a fifteen 
minute cloud passage. As the cloud starts its passage, the solar input power 
goes to zero, the receiver lid closes and the subsystem temperatures start to 
decrease. The helium working fluid pressure (P) will decrease until the pre- 
selected minimum value is reached resulting in engine operation at a low or 
"idle" power level. During the fifteen minutes of cloud passage the receiver 
pool and engine head temperatures decrease while the engine "idle" power 
remains relatively constant. When the cloud has passed, the solar power 
returns and the receiver lid is opened allowing the pool and head temperatures 
to increase. Approximately four minutes is required to heat the subsystem 
back up to the steady-state operating temperature and return to the "equilibrium" 
power level. The reheat time in this case is one minute less than for the 
previously presented normal start-up time because the solar input power is 
greater for this mid-day case. Thus, a fifteen minute cloud passage causes 
the baseline system to operate at "idle" (and minimum "on-line") conditions 
for only nineteen minutes before normal conditions are resumed. 
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• 30 MINUTE CLOUD PASSAGE - BASELINE SUBSYSTEM 

Figure 12 presents the temperature response of the baseline subsystem for a 
thirty minute cloud passage. This case is identical to the previous case 
for the first fifteen minutes of cloud passage. After this period, the head 
temperature reaches the lower operational limit, and engine and valve are 
shut off. The pool and head temperatures slowly decrease in temperature until 
the sun returns. When the cloud has passed (and the solar input power re- 
turns) the receiver lid and valve are opened. The pool and head temperatures 
quickly increase to the engine start level (less than one-half minute) where 
the engine operates at "idle" conditions. From this point on the thirty 
minute case is identical to the previously described fifteen minute cloud 
passage case. 

• 15 MINUTE CLOUD PASSAGE - BASELINE SUBSYSTEM WITH 100 KG OF BUFFER STORAGE 

The temperature response to a fifteen minute cloud passage for the subsystem 
with 100 kg of buffer storage is presented in Figure 13. The system will 
operate at "rated" conditions for approximately six minutes from the latent 
heat of fusion of the eutectic salt. After the buffer storage has completed 
its constant temperature phase change, the engine will operate at "idle" 
conditions. The steady-state operating temperature is reached approximately 
four minutes after the solar power returns. The solar input power that enters 
the receiver which is in excess of that required to operate the Stirling engine 
at "rated" conditions is used to melt the eutectic salt storage. Thirteen and 
one-half minutes of "rated" power operation is required before the engine re- 
sumes operating at the "equilibrium" power level. Thus, for a fifteen minute 
cloud passage, approximately thirteen minutes "idle" opeation occurs for the 
subsystem with 100 kg of buffer storage. 

• 15 MINUTE CLOUD PASSAGE - BASELINE SUBSYSTEM WITH 250 KG OF BUFFER STORAGE 

Figure 14 presents the temperature response to a fifteen minuta cloud passage 
for a subsystem containing 250 kg of buffer storage. For this case, the 
engine will operate at "rated" power for the entire fifteen minutes of cloud 
cover. However, an additional 33.5 minutes of "rated" power operation is 
required, once the solar power continues, while the eutectic salt storage 
"recharges" (melts). 
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• BUFFER STORAGE 

Eutectic salt buffer storage within the receiver subsystem can effectively 
maintain efficient engine operation during moderate solar drop-out periods. 

This is accomplished when the salt gives up its stored thermal energy during 
a constant temperature liquid to solid phase change. A salt with a melting 
temperature near the steady-state operating temperature is desirable since 
the Stirling engine operates less efficiently at lower head temperatures. 
(Reference 4) The analyses of References 7, 8 and 9 suggest NaF/MgF,, eutectic 
salt as being an excellent candidate for thermal storage. The selected mole 
percentages of NaF and MgF^ are determined by the desired eutectic melting 
point. A 75NaF/25MgF ? eutectic salt has been used for the present analyses. 

The effects of adding thermal buffer storage are presented in Figures 15 and 
16. These two figures are based upon the analyses of Reference 7 and updated 
to reflect the effects of operating the Stirling engine at 1800 rpm. 

Figure 15 presents the sodium temperature response during engine operation 
at "rated" and "idle" conditions. Operation was for a subsystem with and 
without buffer storage as well as for an open and closed receiver lid. The 
effect of closing the lid is negligible during "rated" power operation since 
heat lost by radiation and convection out the receiver aperture is small 
compared to the energy removed by the engine. However, with the engine oper- 
ating at "idle" power, the reradiation and convection losses are very signifi- 
cant. For example, by simply closing the receiver lid during "idle" conditions, 
engine operation at 800°C with 100 kg of buffer storage can be extended from 
33 minutes to 60 minutes. 


Figure 16 presents the engine operating time with a closed receiver lid as 
a function of buffer storage mass. As shown, engine operating time during 
the salt phase change can be increased by 6.0 minutes for each 100 kg of 
75 NaF/25MgF2 eutectic salt added to the subsystem. An additional 7.0 minutes 
of engine operation can be obtained from the entire subsystem sensible heat. 
However, it may not be desirable to operate at "rated" power below the steady- 
state head temperature. This is due to the fact that engine efficiency will 
decrease and warm-up time will increase as the engine head temperature decreases. 



Engine operating time is thus gained at the expense of subsystem weight. 

For example, 250 kg of eutectic salt is required to maintain "rated" output 
power for a 15 minute solar drop-out period. Although not considered in 
the present analysis, the weight of encapsulation must also be added to that 
of the salt. It is desirable to use thin wall containers to minimize structural 
weight, cost and temperature drops to and from the salt. Fortunately, by in- 
corporating the following two techniques, additional wall thickness for cor- 
rosion allowance and pressure differentials is not required. 1) By purifying 
the salt and using a "getter",e.g. aluminum, the corrosion of 316 stainless 
steel can be reduced to nil (References 6, 10 and 11). 2) By adding a small 

amount of sodium to the evacuated salt container, the sodium vapor pressure 
inside and outside the container are (automatically) nearly equal all the 
time. (Reference 6) 

"On-line" operating time at "idle" power can be extended by two other methods. 
The first method is by simply increasing the mass of sodium within the re- 
ceiver. The second is by reducing the "idle" shaft output power level from 
the present 10 kW level to the minimum level of approximately 1 kW. This 
would increase "idle" operating time by a factor of three or four. 

However, extending the engine operating time beyond that which is possible 
with the inherent thermal capacitance may not be necessary. A cursory in- 
vestigation of both the Barstow and Lancaster insolation data indicate that 
short cloud passages (less than 3 minutes) occurred more than a thousand times 
at these locations during 1976. Cloud passages on the order of 15 and 30 
minutes in length occurred much less frequently (less than 100 times). There 
were eighteen days in which the total available energy was between 50 percent 
and 10 percent of the annual mean level. Only eleven days had a total energy 
level, less than 10 percent of the annual mean level. Although it is very 
difficult to determine the frequency and length of cloud cover from 15 minute 
insolation data, the above indicates that the duration of a "typical" cloud 
passage is less than 15 minutes. Thus, buffer storage for a longer period of 
time may be unnecessary. 
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• EMERGENCY OPERATION 


Figure 17 presents the receiver wall temperature response during an emergency 
operation in which the engine load suddenly goes to zero. For the conditions 
listed, the wall temperature increases about 65°C per minute for the baseline 
subsystem and about 25°C per minute when 100 kg of buffer storage is added to 
the subsystem. These temperature rates are sufficiently slow to allow de- 
focusing before damage to the receiver can occur. 

• SUBSYSTEM RESPONSE CHARACTERISTICS 

Figure 18 presents response characteristics for the baseline subsystem during 
start-up, in addition to those previously presented in Figure 10. In the 
morning, prior to focusing upon the sun, the receiver lid and vapor valve are 
opened. The receiver wall and sodium pool temperatures rapidly increase when 
the solar input power first comes "on”. The receiver wall temperature is an 
average of 23°C higher than the pool temperature during the fifteen minute 
frozen pool start-up period. The pool and head pressures correspond to the 
saturated sodium vapor pressure at the respective temperatures. As shown, 
the pressure magnitudes and differentials are negligible for temperatures 
below approximately 400°C. As the sodium pool temperature increases above 
this temperature, the pressure differential and vapor density become great 
enough to significantly increase the vapor flow to the engine heater head. 

As a result, the head temperature begins to increase, and the slope, of the 
pool temperature curve decreases. This slope change is a result of the incoming 
solar power now heating both the pool and head thermal capacitances. The 
temperatures, pressures, and vapor mass flow rate increase as the solar power 
continues to enter the system. The Stirling engine is started and operated 
at "idle" power when the head temperature reaches the start temperature of 
450°C. An inflection of the head temperature occurs at this point since some 
of the power which would normally heat the head is now used to drive the 
engine. The sodium vapor mass flow rate peaks, and then decreases as the 
power required by the engine heater head (to increase temperature and drive 
the engine at "idle" conditions) is reduced. The pool and head temperatures 
and vapor pressures continue to increase until the 800°C steady-state operating 
temperature is reached. At this point the engine is controlled to maintain 
constant head temperature. The helium working pressure is increased, and the 
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engine is operated at the "equilibrium" shaft output power level. The power 
that was used to increase the pool and head temperatures is now used for input 
to the Stirling engine. Thus, the sodium vapor mass flow rate also increases 
to the "equilibrium" conditions. 


4.0 CONCLUSIONS 

The development of an improved transient thermal model for the baseline 
receiver-thermal transport subsystem operating with the USS P-75 Stirling 
engine has been presented. The two governing differential equations were 
solved numerically, and detailed results for various transient cases were 
presented. These results demonstrate that the transient performance of the 
subsystem is stable and well behaved, and has long time constants. Salient 
results of the analyses are as follows: 

1) For normal morning start-up, only eight minutes are required to heat 
the system up to the 800°C steady-state operating temperature. In 
the first three minutes the ambient head temperature is raised up 

to the 450°C start temperature. 

2) Normal start-up to "idle" power requires approximately four minutes 
for the subsystem with 100 kg of buffer storage. The additional time 
required to reach "equilibrium" power (and to melt all of the salt) 
is a function of the salt melt temperature and mass. 

3) For start-up with a frozen pool, 10 and 15 minutes are required to 
reach the start and steady-state operating temperatures, respectively. 
Less than one minute is required to melt the pool. 

4) Due to the inherent thermal capacitance of the subsystem, the engine 
will operate at "idle" power for fifteen minutes after the solar power 
goes to zero (evening or cloud passage) . 
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5) If the receiver lid remains closed during extended "bad" weather con- 
ditions, the pool temperature is predicted to stay above freezing in 
excess of three days. With the receiver lid open the pool temperature 
drops to the sodium freeze temperature in about sixteen hours. 

6) The baseline system will stay "on-line" during a fifteen minute cloud 
passage without adding buffer storage. 

7) Approximately 250 kg of buffer storage is required to maintain "rated" 
power throughout a fifteen minute cloud passage. 

8) During an emergency operation, the temperature rise rate is sufficiently 
slow to allow defocusing before damage to the receiver can occur. 

The current transient thermal model is complete, except for any minor changes 
that may be required as a result of improvements to the receiver, thermal 
transport duct, and Stirling engine. Additional transient cases will be run 
for any changes to the baseline system. 



TABLE I 


Receiver-Thermal Transport Subsystem Data 

Steady state operating temperature = 800°C (1472°F) 

Engine start temperature * 450°C (842°F) 

■ 15 ° c <59 ° F) 

- -10°C ( 14 °F) 

W-N. * ,8 ° C < 208 ° F > 

LI-SALT ” 800 ° C <U72 ° F) 

Sodium mass = 33 kg (73 lbm) 

Structural mass = 150 kg (331 lbm) 

Effective head mass = 35 kg (77 lbm) 

Salt mass * variable 

Salt heat of fusion - 10.578 kW - min/kg (273 Btu/lbm) 

Duct diameter - 0.0762m (3.0 inch) 

Total transport loss coefficient =7.0 

Insulation conductivity = .03634 W/m-°K (.021 Btu/hr-ft-°F) 

- 6.389 x 10" 12 kW/°K 4 (2.077 x 10' 9 Btu/hr-°R 4 ) 

- 6.90 x 10 -4 kW/°K (1.308 Btu/hr-°F) 

P MIN = 2.5 MPa (25 atm) 

P VAV = 15.0 MPa (150 atm) 
maa 
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FIGURE 4 
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FIGURE 5 

SODIUM GAS CONSTANT AND COMPRESSIBILITY FACTOR 
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FIGURE 6 

THERMAL RESPONSE FOR NORMAL START UP FROM PREVIOUS DAY 
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FIGURE 8 

THERMAL RESPONSE FOR NORMAL START UP FROM PREVIOUS DAY WITH BUFFER STORAGE 
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FIGURE 9 


THERMAL RESPONSE FOR NIGHTTIME AND EXTENDED "BAD" WEATHER 
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FIGURE 10 

THERMAL RESPONSE FOR START UP WITH FROZEN POOL 






FIGURE 12 


SUBSYSTEM TEMPERATURE RESPONSE TO A 30-MINUTE CLOUD PASSAGE 


N * 1800 RPM 
NO BUFFER STORAGE 


CLOUD PASSAGE 


SOIAR 


SOLAR 


"EQUILIBRIUM 1 

ENGINE SHAFT 
OUTPUT POWER 




11:30 11 


LID 

CLOSED 

VALVE CLOSED 
ENGINE OFF 

LID & VALVE 
OPENED 

ENGINE START UP 

SOIAR TIME (HR:MIN) 








ENGINE SHAFT 
OUTPUT POWER 









FIGURE 15 


SODIUM TEMPERATURE RESPONSE AT RATED & IDUi ENGINE POWER 
N-1800 RFM 

S^-Engine at rated power, no buffer storage S.^-Engine at Idle power, no buffer storage 

S^-Englne at rated power, 100 Kg buffer storage S^*Englne at Idle power, 100 Kg buffer storage 






P -65KW SHAFT 


P-150 AIM 


p 2 -ikw SHAFT 


33 Kg SODIUM 

150 Kg STAINLESS STEEL 

SALT-75 NaF/25MgF 2 


LID CLOSED 

LID OPEN 


ENGINE OPERATING TIME (MIN) 







FIGURE 16 

ENGINE OPERATING TINE Willi BUFFER STORAGE 
N-1800 RPM 
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FIGURE 18 

SUBSYSTEM RESPONSE DURING START-UP 
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APPENDIX I: NOMENCLATURE AND UNITS 


‘SlSH 
SUCT 
C SALI 
C HEAD 
C SODIUM 
C STRUCTURE 

°» 

f abs 

r DROP 


f ENG 


f 


LID 


^ALVE 


g 

g cool 

g eng 


2 

collector area (m )* 

2 

duct cross-sectional area (m“)* 

buffer storage capacitar.ce (kW - min/K)* 

engine head thermal capacitance (kW - min/K)* 

sodium thermal capacitance (kw - min/K)* 

pool and duct structure thermal capacitance (kW - min/K)* 

sodium vapor characteristic velocity (m/sec) 

clear sky absorption path length correction factor 

solar dropout factor* 

1 0 for cloud passage or night 
1 for clear sky 

engine factor* 

! 0 for engine off 
1 for engine running 

receiver lid factor* 

! 0 for lid closed 
1 for lid open 

sodium valve factor* 

! 0 for valve closed 
1 for valve open 

1 kg m/sec 2 /N 

effective heat transfer coefficient for engine cooling 
system (kW/K)* 

conductance for heat loss through the engine when the 
engine is off (kW/K)* 
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u 


HEAD 

G L1D 

"Xa 

g pool 

g wall 

?**> 


K 

L 

l 

M 

m 


ABS 


m v 

<v 


p 

P HEAD 

P POOL 

P REF 

P v 

(PR) 


AMB 


(PR) 


HEAD 


Q CONDH 

^CONDP 

^CONV 

^ENG 


* conductance for engine head conduction loss (kW/K)* 

- lid conductance in closed position (kW/K)* 

* sodium pool boiling convection coefficient (kw/K)* 

=■ conductance for receiver conduction loss (kw/K)* 

■ receiver cavity wall conductance (kW/K)* 

=• convection loss coefficient (kW/K)* 

2 

* solar flux above earth's atmosphere (kW/m“)* 

■ clear sky absorption for 90° solar elevation* 

* site latitude* 

* pressure exponent for engine input power* 

■ sodium vapor mean molecular weight (kg/kg-mole)* 

* vapor pressure exponent* 

■ sodium vapor flow rate (kg/min)* 

- sodium vapor flow rate for choked flow through the duct 
(kg/min) 

- engine mean pressure level (atm)* 

* sodium vapor pressure at engine head (Pa) 

-* sodium vapor pressure at the pool (Pa) 

* reference vapor pressure (Pa)* 

» sodium vapor pressure (Pa) 

* input power ratio for coolant temperature different 
from 15°C 

- input power ratio for head temperature different 
from 800°C 

■ head external conduction loss (kW) 

» receiver and duct external conduction loss (kW) 

- receiver internal convection loss (kW) 

■ power input to the engine (kW) 


L-39 


TffiAD 

w 
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^OUT 

Q POOL 

q rr 

^SOLAR 

^TRAN'S 

Q 100 

R 

t 

t amb 

t cool 

t head 

t pool 

t ref 

T 

SKY 

T 

v 

T 

w 

Z 

3 

Y 

r 
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net head power (kw) 

reference engine output power at zero engine pressure 

(kW)* 

engine output power (kW) 

net pool power (kW) 

receiver re-radiation loss (kW) 

solar heating rate at the receiver (kW) 

power transmitted through the vapor duct (kW) 

reference engine input power at 800°C head temperature 
and 100 atm engine pressure (kW)* 

universal gas constant 

8314.34 (W-sec/kg-mole-K) 

time variable (min, solar) 

ambient air temperature (K)* 

mean temperature of engine coolant (K) 

engine head temperature (K) 

sodium pool bulk temperature (K) 

reference temperature for vapor pressure function (K)* 
sink temperature for re-radiation (K)* 
vapor pressure constant (K)* 
receiver cavity surface temperature (K) 

= vapor compressibility factor* 

* slope of engine output vs. input power at 800°C head temp. 
■ sodium vapor ratio of specific heats* 
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DSH 

'ENG 

<"V 



-* solar declination angle* 

* sodium heat of vaporization (kW - min/kg)* 
■ collector efficiency* 

= engine efficiency 

a 

= re-radiation coefficient (kW/K )* 



denotes input 






APPENDIX M 

SHADOWING ANALYSIS FOR PARABOLIC 
CONCENTRATOR FIELDS 

1.0 INTRODUCTION 

This Appendix presents the development and results of a shading/shadowing* 
analysis for a field of parabolic concentrators. A loss of energy occurs 
whenever solar insolation is prevented from reaching the aperture of a 
collector. This solar blockage can occur as a result of one (or more) 
collectors shading another. The percentage of the annual energy which is 
lost is a function of both the amount of energy available and the percent 
shading at any given instant which in turn are functions of the site 
s location and the solar time. The amount of shading is also related to the 

field layout and the collector spacing. Thus, the objective of this study 
is to evaluate the annual energy loss for various field layouts and collector 
spacings. The results presented in this study are based upon the 15 minute 
solar insolation data for 3arstow, California for the entire 1976 year. 

These results were used to optimize the collector arrangement based upon 
site location, land costs and costs dependent upon plant layout (for example, 
electrical cabling, fencing, etc.) 

2.0 METHODS OF ANALYSIS 

The collector shading, as a percentage of the total collector area, is 
derived from the relationships presented in Section 4. The percent shading 
is determined for any field of parabolic concentrators given the following 
information: 

NNS « Number of concentrators in the North-South direction. 

NEW * Number of concentrators in the East-West direction. 

LNS =* Ratio of the concentrator spacing (center- to-center) to concen- 
trator diameter in the North-South direction. 

LEW ■ Ratio of the concentrator spacing (center-to-center) to concen- 
trator diameter in the East-West direction. 

9 * Aximuth angle of the sun (referenced from North) . 

A 

cp E ■ Elevation angle of the sun (referenced from the horizon) . 

The values of cp £ and co were calculated for each 15 minute time period using 
the equations presented in Section 5.0. 


^Shading and shadowing are used interchangibly in the Appendix. 
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The percent energy loss for fixed values of field layout and sun angles is 


,OSS AVG CONC d,n 


inhere 


Average shaded area for the entire concentrator field 


Concentrator area (unity) 


CONC 


Solar insolation (direct, normal) 


The values of the percent energy loss were compared with the percent loss in 
Capacity Factor for the baseline system using the P-75 Stirling engine. The 
difference between the two methods was sufficiently small (less than 3% of 
the respective values) to permit the use of this approximation. The percent 
energy loss was calculated for greater than or equal to 10 . This corres- 
ponds to the start-up and shutdown elevation of the baseline system. If the 
start-up and shutdown elevation was zero, the percent annual energy loss 
would increase by only about two percent. This is because the solar insolation 
during the first and last hour of a day are quite low. 


The collector packing fraction was determined as follows 


CONC 


Area of the solar field 
( (NNS-1) *LNS+1) * ( (NEW- 1) *LEW4 1) 


3.0 RESULTS 


Figures M-l through M-4 present the percent annual energy loss due to collector 
shadowing for various plant layouts of an 18 collector system. 


Figure M-l presents the percent annual energy loss from collector shading, 
as a function of packing fraction and collector spacing, for a 6 by 3 plant 
layout. As expected, the energy loss due to shading decreases as the packing 
fraction decreases (collector spacing increases). However, for a given 
packing fraction and field configuration there is an optimum collector spacing 
Increasing the North- South collector spacing (LNS) beyond about 2 dish 
diameters is not required since the shadowing (and hence energy loss) in this 
direction is negligible at greater distances. Likewise, a 5 dish diameter 
spacir.g in the East-West direction (LEW) is sufficient to reduce the shadow in 
this direction to a negligible level. Various combinations of LNS and LEW 
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spacing produce the minimum energy loss at a given packing fraction. For 
example, the baseline system (6x3 matrix) with a packing fraction of about 
25 percent has a minimum annual energy loss of less than 1 percent. This 
occurs for collector spacings of approximately 1.5 and 3.0 dish diameters in 
the North-South and East-West directions, respectively. Similar results are 
obtained for other collector field arrangements as shown in Figures M-2 
through M-4. 


4.0 DERIVATION OF THE SHADING RELATIONSHIPS 


This section presents the method to determine the percentage of shading that 
a solar collector experiences as a function of collector spacing and sun 
position. Figure M-5, which shows the top and side views of a cluster of 
four collectors, defines the geometry of this analysis. For certain combina- 
tions of sun azimuth angle (<? ) , sun elevation angle (<?„) , and collector 
spacing (L™, and L„ ) , collector 3 can experience different types of shading 
from neighboring collectors. Some combinations may yield individual shading 
by collector 4 (called East-West or 4-3 shading) , collector 2 (called North- 
South or 2-3 shading) , or collector 1 (called diagonal or 1-3 shading) . 

Other combinations may produce cross shading when collector 3 is shaded 
collectively by collectors 1 and 2 (called 1, 2-3 shading) or collectors 
1 and 4 (called 1, 4-3 shading). This analysis examines each of these 
possibilities, and determines the extent of shading that each produces. 

a. East-West Shading . Consider first the East-West shading (for example, 
the shading of collector 4 onto collector 3) as illustrated in Figure M-6. 

*4-3 ' hv ,ln( V 

' D/ EW 

*4-3 ' COa ®A> 

- EH cos <V 

where 

^EW " 

D is the dish diameter 

See Figure M-7 

*’ 4 -3 ’ *4-3 3iD 

- D l EW sin (0 A ) sin (0 £ ) 

See Figure M-8 

d 4-3 A/*'4.3 > 2 * <W* 

- d/j^ / sin 2 (0 A ) Sin 2 (0 £ ) + cos 2 «* A ) 
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The common area, A, , , is given by: 

2 

A 4-3 * T (6 1 “ 3in •l 5 

where the angle, 9,, con be found as follows: 
*L ' 

2rcos t'— - d, , 

L 4-3 


( 1 ) 


9^-2 cos 


-1/ 4-3 


C) 


The following test is used to determine whether or not 4-3 (or East-West) 
shading exists. 


If 


d 4-3 


>1; no shading occurs 

A 4-3 ’ °» d i * 0 

<1 ; shading occurs 

- eq.(l); - eq. (2) 


b. North- South Shading . Figure M-9 illustrates the case of North-South 
(or 2-3) shading. The formulation of shading for this case is similar to that 
used in the East-We3t shading case. 

*2-3 ' 'US coa <V 

- D *NS cos <V 

^2-3 ■ Ss «A> 

• 0/ ss sln «A> 


where 


^ NS " SlS^ 


See Figure M-10 


* 2-3 " *2-3 Sln 
See Figure M-ll 


- D^ ns cos (0^) sin (0 E ) 
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(y 2-3 ) 


P/ NS\ COs2 Sin2 + Sin2 


The common area, ^2-3* * ,a Siven by; 

A 2 _ 3 - D 2 (0 2 - sin 9 2 ) 
where the angle, 9,,, is given by; 


9 0 - 2 cos 


-1/d 


2-3 


(3) 


(4) 
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The following test is applied to determine whether or not 2-3 (or North-South) 
shading occurs. 


>1; no shading occurs 


shading occurs 

A, , * eq.(3); * eq.(4) 


c. Diagonal Shading . Consider now the diagonal shading (for example 
shading of collector 1 onto collector 3) as shown in Figure M-12. 





The common area, A^_ 3> is given by; 

A 1 _ 3 = (© 3 - sin 63 ) 


where 


2 cos 


The following test is used to determine whether or not diagonal (or 1-3) 
shading exists. 


> 1 ; no shading occurs 
. A. » 0; 0 ■ 0 

If till 1-3 3 

< 1 ; shading occurs 

Al _3 - eq (5) ; 0 3 - eq (6) 

d. Cross Shading; . Cross shading on any given collector occurs whenever 
two shadows, cast by other collectors, intersect and share a common area. 
Figure >1-15 shows how collectors 1 and 2 may cross shade collector 3 and how 
collectors 1 and 4 may cross shade collector 3 . 

Figure M-16 is an enlarged drawing of the condition where collector 3 is 
double shaded by collectors 1 and 2 . 


From the law of cosines; 


cos” 1 / d 2 + d 2 - d 2 

l 1-2 1-3 2-3 

V 2 d l -2 d l-3 ' 

\ 2 d 2-3 d l-2 


- 1 / 222 
S 3 - C ° S f d 1-3 + d 2-3 - d 1-2 

\ 2 d l-3 d 2-3 / 


where 


From the geometry of the figure, it can be shown that, 

- % (® 3 + e 2 ) - 6 l 
y 2 - h <e 2 + e L ) - 5 2 
y 3 - % (^ + « 3 ) - 63 
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Ihe areas of the segments S^, S^, and are given by; 

Sj. * D_2 (Y - sin Y^) 

8 

S 5 * £1 (Y, “ sin Y») 

8 4 

S T “ 2l (Y, - sin Y,) 

8 

The chord lengths b^, b 2 , and b^ are expressed by, 
b L » D sin (%y l ) 

b 2 =• D sin (Yy^) 

b 3 =■ D sin (%y 3 ) 

The area of the triangle enclosed by chords b^, b 2> and b 3 is; 

A tri - yj. x (x - b L ) ( x - b 2 ) (x - b 3 ) 

where 

x =» h (b L + b 2 + b 3 ) 

The area of the region common to all three circles, 2 3 is given as; 

*1, 2-3 ' S i + S 2 + S 3 + Vi < 7 > 

The following test is used to determine the extent of cross shading onto col- 
lector 3 due to collectors 1 and 2. 

< 0; no cross shading exists 

" V 2-3 ■ 0 

> 0; cross shading exists 

* 1 , 2-3 ■ ec >- (7> 

The analysis for 1, 4-3 diagonal shading is identical to that of 1, 2-4 shading 
except d^_^ is used in place of d^_ 2 , and d^_ 3 is used in place of d 2 3> 

5.0 SOLAR ANGLE RELATIONSHIPS 

The azimuth and elevation angles of the sun are determined by the following 
expressions*: 

*Anon., "On The Nature and Distribution of Solar Radiation", Prepared for the 
U.‘S.* Department of Energy, HCP/T2552-01, March 1978. 
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LcosL cos6 cosh + sinu sinoj 


r ) 


where 


' u 

r A 


sin 


* cos" 1- _ (sins - sinL si: OE )y/(cosL 
h - 15 (H g - 12) 

5 - 23.45 sin L1.008(n-80)] 1 
5 * 23.45 sin L .965(n-80)] 81 
5 * -23.45 sin l .975(n-266)J 267 


cos©,) 


< n < 80 

< n < 266 

< n < 366 


L * site latitude (degrees) 

S « solar declination (degrees) 
h * hour angle about noon (degrees) 
H s » true solar time (hours) 
n ■ number of the day of the year 
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1.0 INTRODUCTION AND BACKGROUND 

w 

The baseline receiver coupled with the Stirling Cycle Engine utilizes 
metallic sodium (Na) as a heat- transfer medium. The excellent physical pro- 
perties of the liquid-metal (sodium) make it possible to design a simple and 
efficient heat-exchanger (boiler) which transfers concentrated sun energy into 
hot sodium vapor. This hot vapor at less than atmospheric pressure is then 
used to operate a Stirling Engine. The high heat transfer coefficient of 
sodium and its relatively low operating pressure at high temperatures are 
..nly two of the favorable properties. 

Figure 1 shows a sketch of the baseline system. The major components of 
the system are receiver core filled with sodium, cavity door, butterfly blocking 
valve, expansion bellows, vapor pipe, Stirling Engine, and generator. The 
entire system up to and including the engine head is covered with thermal in- 
sulation. The receiver core contains approximately 33 Kg (73 pounds) of sodium 
which is heated to approximately 800°C. The hot sodium vapor from the receiver 
core is ducted by a vapor pipe through a butterfly valve and bellows to the 
engine head. During the night-time shut-down periods the butterfly valve is 
closed to block vapor flow to the engine, and the cavity door is closed to 
minimize heat loss. The complete baseline system shown in Figure 1 has a total 
weight of approximately 1088.6 Kg (2400 pounds). 

This report presents summarized results from an in-depth study of metallic 
sodium. The topics presented are commercial applications of metallic sodium, 
hardware found in sodium systems, corrosion of materials (stainless steel and 
nickel-base alloys) exposed to liquid sodium, and sodium handling. 

2.0 COMMERCIAL APPLICATIONS OF METALLIC SODIUM 

Presently, the applications of metallic sodium are divided into two major 
categories — those based on the physical properties of sodium and those based 
on the chemical reactions of sodium. This report will be confined to the uses 
based on physical properties and will neglect the uses based on inorganic and 
organic chemical reactions. 

Probably the best known use of sodium is as a heat transfer medium. This 
application is very evident if the heat transfer coefficients (boiling, condens- 
ing, and convection) of sodium are compared to other working fluids. The ability 
of the sodium to rapidly transfer heat at high temperatures makes it very useful 
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for solar heating applications where high temperature transients are encountered. 
Other unique properties of sodium which make it a very desirable working fluid 
includes: 

a) Low operating vapor pressure over a wide temperature range. (This 
quality is very important for mechanical design and fabrication of 
heat exchangers) . 

b) Sodium has previously been used in a Stirling engine. (Reference 1) 

c) High boiling point. 

d) Relatively high merit number over a wide temperature range. 

A few examples follow which describe commercial applications of metallic 
sodium. 

2.1 ENGINE VALVES 

For many years sodium has been used as a heat transfer medium in sodium- 
cooled valves in trucks, cars, and aircraft engines. Sodium is contained within 
a hollow space inside the valve. Heat is transferred from the head of the valve 
to the stem by the splashing of the liquid sodium as the valve moves up and down 
during engine operation. Figure 2 shows operating temperatures for a conventional 
solid stem valve and a hollow stem sodium- filled valve. Figure 2 shows that with 
60 percent of the hollow valve stem filled with sodium the operating temperature 
of the valve neck and stem is reduced approximately 190°C. This substantial 
reduction of operating temperature of the valve material prolongs valve life. 

2.2 NUCLEAR POWER PLANTS 

Because of the very desirable properties of sodium, it is an ideal heat 
transfer medium for nuclear power plants. These properties include: 

• Low neutron absorption 

• High boiling point 

• High heat-transfer efficiency 

• Extreme thermal stability 

• Ease of pumping 

• Freedom from corrosion in the absence of impurities 

• High temperature operation at low vapor pressure 
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Typically, the heat removed from the reactor core is transferred to a 
steam generator through two separate sodium loops; a primary and secondary system. 
This dual system is necessary to prevent radioactive sodium in the primary loop 
from coming into contact with the fluid used to drive the steam generator. There- 
fore, the radioactive primary loop is separated from the water system by means 
of a secondary sodium system. 

For several years a sodium-cooled nuclear reactor built by Atomics Interna- 
tional at Santa Susanna, California, has supplied heat to generate electricity 
for a commercial grid. Because liquid metal (sodium) has such good heat-transfer 
capability and low operating pressures at high temperature, this nuclear power 
plant is operated safely and efficiently. 

A few of the many other non-chemical applications of liquid sodium are 
listed below: 

• Liquid sodium has been used as a heat-treating medium to provide uniform 
temperature distributions. 

• In the die-casting of magnesium, sodium is used for core cooling of 
casting machines. 

• Another instance of the successful application of sodium to an industrial 
heat-transfer problem is in shale oil retorting. Special plows which 
distribute shale are exposed to 1200°C operating temperatures and cooled 
by sodium. 

• For the control of glass mold temperatures. Sodium is used to absorb 
heat during the glass molding cycle, and also preheat empty molds. 

• Liquid metal (sodium) is used in the reheat cycle of large steam plants. 
Liquid metal heaters, located in the boiler, carry hot sodium to the 
second stage turbine reheater. Superheated steam, which has expanded 
through a high pressure primary stage turbine, is then reheated before 
passage through the second turbine. 


3.0 SODIUM SYSTEMS HARDWARE 


3.1 PUMPING /TRANSFERRING LIQUID SODIUM 

Both mechanical and electromagnetic pumps have been used with relatively 
good success as a means of pumping molten sodium. The primary problem in 
designing pumps used to trsmsfer sodium (as with valves) is to provide satis- 
factory seals. Currently there are various types of pump seals being used - 




a frozen seal, gas* sealed, and mechanical seals. Examples of these types of 
seals are shown in Figure 3. 

Electromagnetic pumps have been used successfully for liquid sodium recir- 
culation systems. These pumps are totally sealed with no moving parts, and 
provide long periods of trouble-free service. A wide range of pump models 
providing different flow rates and pressure outputs are available. Figure 4 
shows a diagram of a D-C Faraday-Type pump. 

Mechanical pumps that have been used successfully to transfer liquid 
sodium are the gas-sealed submerged centrifugal pumps, piston and diaphram pumps, 
and rotary gear-type pumps. These pumps have high operating efficiencies and 
are capable of maintaining high system pressure for continuous service. A 
large number of pumps with a variety of flow rates and pressure outputs are 
available from commercial pump companies. A single stage centrifugal pump and 
a gas-sealed type centrifugal pump are shown in Figure 5. 

For many operations, liquid sodium can be transferred by differential- 
pressure methods. One method used for inducing sodium flow is to establish a 
vacuum in the receiving tank or equipment. A low-pressure nitrogen atmosphere 
is maintained in the tank from which the sodium is being withdrawn. The 
pressure differential thus creates sodium flow. A simple air-cooled, recipro- 
cating vacuum pump can be used to produce the necessary vacuum. To prevent 
drawing liquid sodium into the vacuum pump, it is common practice to install 
a barometric loop ahead of the pump. 

3.2 FILTERS/COLD TRAPS 

Raw sodium from the manufacturer contains small quantities of impurities 
such as oxides, chlorides, and calcium metals. However, because the solubility 
of these impurities is reduced at low temperatures, the sodium can be purified 
by filtration, a number of devices and techniques have been employed for removal 
of these impurities from molten sodium. These include: 

• Chemical "getting" of oxygen by calcium or titanium. 

• Filtration. 

• Cold-trapping. 

Chemical purification by "getting agents" is used primarily for oxygen 
content reduction. Two of the more common gettering agents are calcium and 
titanium. The calcium getter material is soluble in sodium and forms an oxide 
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with limited solubility. These oxides are then reacted with impurity materials 
soluble in sodium and removed by filtration. 

The second gettering agent, titanium, Is insoluble in sodium. The titanium 
getter reduces sodium oxides to insoluble getterer oxide which is removed from 
the system by filtration. 

The second technique commonly used to purify sodium is by filtering at 
approximately 150°C through Micro Metallic filters. These filters have a pore 
size of approximately 15-microns and are constructed of sintered stainless 
steel. Thij type of filtering can provide an oxide concent as low as 10 parts 
per million (ppm) . This oxide level was found to be satisfactory for the 
proposed FACC baseline sodium receiver. 

A cold trap is a basic purification device for flowing sodium systems. In 
this unit, sodium is cooled down to a temperature of approximately 170°C. At 
this reduced temperature, oxides and other impurities become insoluble and pre- 
cipitate out of the sodium. Under these conditions, the trapped impurities can 
be permanently removed from the sodium system. In this form of purification, 
the oxide content of a sodium system can be maintained as low as about 5 to 10 
parts per million. This level is satisfactory for most systems constructed of 
stainless steel material. 

A typical cold trap, shown in Figure 6, would consist of a canister con- 
taining 10-15 micron pore 3ize stainless steel wool sized such that the flowing 
sodium would be kept at approximately 170°C. Maintenance of a cold trap would 
involve periodically replacing the contaminated stainless steel wool. Since the 
FACC baseline sodium receiver operates isothennally, cold trap devices are not 
necessary for sodium filtering. 


3.3 VALVES 

Typical valves found in sodium systems are blocking valves, throttle-type 
valves, and ball-type check valves. These valves are designed to perform spe- 
cific tasks and meet the general requirements to provide safe and efficient 
operation. 

The general requirements for valves for liquid sodium systems are: 

• Zero valve stem leakage. 

All valves must be specifically constructed for zero leakage to the 
surroundings to prevent impurities, such as oxygen, from entering the 
system. This type of leakage prevention requires all-welded construc- 
tions of valves with either bellows seals or freeze seals. 
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• Corrosion Resistance. 



Corrosion resistance materials are used in all valves. Typically, 
the standard materials used are types 304 and 316 stainless steels. 
However, valve seats, discs, and other wearing parts should be con- 
structed with hard facing materials in order to prevent galling. The 
possibility of intergranular corrosion requires that special pre- 
cautions be taken to insure sound metal parts with no porosity. The 
stringent requirements make the use of forged valve bodies very 
desirable. 

Thermal Shock. 

Valves must be designed to resist the pipe stresses resulting from 
thermal cycling and shock. Precautions should be taken in valve 
design to prevent binding due to thermal expansion. 


A typical ball check valve is shown in Figure 7. This valve permits flow 
in only one direction and is operated by the pressure difference across the 
valve. Figure 8 shows a cross-section view of a pis ton -operated blocking valve. 
The FACC baseline receiver will utilize a butterfly valve similar to the valve 
shown in Figure 9. 







4.0 CORROSION OF MATERIALS 

This section of the report will focus on the corrosion behavior of stain- 
less steels (300 series) and nickel base alloys in high temperature sodium. 

The various parameters which affect the rate of mass transfer, rates of material 
loss, and carbonization of stainless steels, have been studied. These para- 
meters include the effect of temperature in the range 600 to 900°C, the effect 
of impurity content of the liquid metal (oxygen, sulfur, hydrogen, carbon) and 
the effect of alloying content such as high percentages or nickel. The effect 
of fluid flow rate will not be considered because it is not applicable to the 
FACC baseline sodium receiver. 

Each of the above parameters will be discussed separately and comparisons, 
where appropriate, are made between various references. 

4.1 MASS TRANSFER 

A common mechanism of corrosion in liquid-metal systems is mass transport. 
Mass transport results from a coexistence of a temperature differential and an 
appreciable thermal coefficient of solubility. Corrosion deposits from the 
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structure materials will tend to dissolve in regions of high chemical activity 
(hot areas), and will tend to deposit from solution in regions of low chemical 
activity (cold areas) . This continued transfer of material in the system can 
accelerate corrosion attack in the areas of high activity. The precipitated 
material will accumulate in cooler regions of the system and may eventually 
cause plugging of flow tubes. At present most sodium pumped systems utilize 
"cold trapping" as a means of continually filtering the flowing sodium. These 
cold traps, which maintain a temperature of approximately 170°C, provide the 
necessary cold area for material deposit. 

Mass transport of impurity substances virtually does not exist in the 
FACC baseline sodium receiver design. The temperature differential within the 
so dium boiler is not sufficient to cause the impurities held in suspension 
to precipitate out of the sodium. The rate of material loss and corrosion of 
receiver material are therefore minimised. 


4.2 EFFECTS OF TEMPERATURE AND OXYGEN CONTENT 

Typically, corrosion experiments involving high-temperature sodium are 
conducted using electromagnetically pumped sodium loops. These loops employ 
flowmeters to monitor fluid velocity through hot-and-cold-leg sections. Heat 
is supplied to the system by wound heaters placed around the hot-leg section. 
The cold-leg consists of a cold trap which collects material deposits during 
the time of system operation. Measurements of mass transfer rates are deter- 
mined by gravimetric and metallographic analysis of the cold trap deposits. 

The corrosion rate data presented in the following paragraphs is based on 
operating systems and measuring techniques described above. 

The mass transfer rates of nickel-base materials (Hastelloy, Inconel) and 
300 series stainless steel are presented in Figure 10 as a function of tempera- 
ture. Figure 10 shows a marked temperature dependence of mass transfer in 
Inconel and Hastelloy materials between 648° and 816°C. The mass transfer rate 
of the nickel-base alloys is very much higher at 816°C compared to 300 series 
stainless steels under comparable conditions. Similar findings were noted by 
Russian investigators (Reference 3) in comparable high temperature evaluations. 

The cold-leg deposits in both the Inconel and Hastelloy systems were 
composed primarily of nickel with small amounts of chromium and virtually no 
iron. The depth of intergranular penetration measured after 2000 hours between 
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704 and 816°C was 1 to 2 mils. Extrapolation of time would result with approx- 
imately 6.57 mils/year material loss. 

The deposits found in the stainless steel systems were composed predom- 
inantly of chromium with small equal amounts of iron and nickel. The steels 
containing 187. chromium and 107. nickel (types 304 and 316) were the most cor- 
rosion resistant of the various types tested. Test results from Reference 2 
indicated a corrosion rate at 816°C of approximately 2.92 mils/year. 

Of equal importance when selecting materials for sodium service concerns 
the effect of oxygen contamination on corrosion. Contrary to the high mass 
transfer rates of nickel-base alloys resulting from high temperature levels 
(816°C), increases in oxygen content levels up to 500 ppm had no measurable 
affect on the material corrosion rates. However, for stainless steel, the 
oxygen content of the sodium is one of the most important parameters. Test 
results (Reference 3) indicate that if the oxygen level is reduced from 25 to 
10 ppm the corrosion rate is reduced by 50% at 650°C. For an oxygen content 
level of 10 ppm the corrosion rate of stainless steel is approximately 2.64 
mils /year. (This material loss rate was used for the FACC baseline sodium 
receiver analysis at an operating temperature of 816°C.) Therefore, this 
study implies that a greater latitude in sodium purity can be permitted for 
nickel-base alloys than for stainless steels. 

Figure 11 presents a comparison of corrosion rates (mils/year) for 300 
series stainless steels from References 1 and 2. For an oxygen content of 
10 ppm and a temperature range from 700 to 816°C, good data correlation can 
be seen. The corrosion rate curves from the two references at 25 ppm oxygen 
content correlate at 650°C only, and differ by 8 mils/year at an operating 
temperature of 816°C. An oxygen content of 10 ppm was assumed for the base- 
line receiver design analysis. 


4.3 CARBURIZATION OF STAINLESS STEEL 

Research studies indicate that a reduction of room temperature ductility 
of 304 and 316 stainless steel is likely with increases in carbon content 
when the materials are exposed to liquid sodium. This carbon increase can 
be attributed to the following: 

• Carbon dioxide (CO 2 ) impurities in the system are reduced by the 
molten sodium to form carbon and oxides. 
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• Higher than normal carbon percentages present in the material following 
manufacture. 

• Small amounts of carbon present in the sodium. 

Long term sodium loop tests (Reference 2) at 650°C show that the average 
carbon levels in both stainless materials (304 and 316) increased approximately 
.25%. Both material specimens contained carbon gradients through their cross- 
sections. The carbon levels that resulted did not materially affect the high 
temperature ductility of the materials. The type of carbide resulting from the 
carburization of the steel surface was identified as Fe^Jk^C, which is dispersed 
in the steel matrix. 

By contrast, data presented in Reference 3 showed only small increases in 
the carbon content. For long term exposure (1300 hours) the exposed samples 
indicated that the carbon content increased slightly to approximately .036%. 
Further testing (Reference 4) was conducted to evaluate the room temperature 
mechanical properties of the stainless steel specimens. These test results 
Indicated that the tensile strength of the exposed sample was above the 
minimum value specified by the ASTM for type 304 stainless steel pipe, whereas 
the yield strength was at the minimum value specified. 

The filling procedure for the baseline sodium receiver is outlined in 
a later section of this report titled "Sodium Handling". The handling, clean- 
ing, and filling procedure described therein will ensure very low levels of 
impurities trapped within the system after filling. Therefore, we are not 
anticipating carburization (decrease in ductility) of the stainless steel 
receiver material. 

5.0 SODIUM HANDLING 

Commercial sodium is available and can be purchased in brick form (1, 25, 

5, 12 pound sizes) and cast solid form. These forms are available in regular 
purity and argon (reactor) purity grades. The argon (reactor) grade sodium 
has fewer impurities than the regular grade, and will be used in the baseline 
receiver tank. 

The brick forms of sodium are packed and delivered in returnable 55-gallon 
steel barrels. These barrels are fitted with airtight covers to maintain a 
protective nitrogen atmosphere. Cast solid sodium is available in 55-gallon 
steel drums and in 5-gallon steel pails. Sodium delivered in these non- 


returnable containers can be heated directly and dispensed as a liquid. 

Figure 12 shows a sodium-tilling system which is currently being used at FACC 
to service experimental sodium/sulfur battery cells. TMs system will be 
utilized to fill the baseline sodium receiver tank. (The filling procedure 
is discussed in the section 5.1 of this report). The closed-filling-system 
illustration shown in Figure 12 reduces handling costs and eliminates oxide 
formation during the filling procedure. The empty 5-gallon sodium pails must 
be cleaned free of sodium before disposal or other use. Tech-Report 015 
describes safety precautions which pertain to container cleaning. 

5.1 RECEIVER PREPARATION AND FILLING 

The corrosion data presented in the previous sections of this report 
indicates the importance of minimizing the oxygen level (ppm) and oxides formed 
when filling the sodium (baseline) receiver. In an effort to minimize these 
corrosion reactants the following surface cleaning treatment and filling pro- 
cedure is recommended: 

1. Vapor degrease component parts of receiver tank prior to welding 
with trichloracthlyene for 45 minutes. 

2. Abrasive cleaning of all welded joints by impinging with Aluminum 
oxide material (#60) . This process should be done during the 
assembly sequence to all welds exposed to sodium. 

3. Nitric acid (technical grade) dip of completely assembled receiver 
tank. 

4. Rinse receiver in distilled water and then in ethylene alcohol. 

5. Dry receiver in oven. 

After the above steps are completed the receiver tank is ready to be 
filled. To eliminate and remove oxides, the following filling procedure is 
recommended : 

(Refer to Figure 12) 

1. Remove trapped air in the system - Close valve D first, then, in 

the following sequence, open valves C, B, and A. This will allow the 
entire system to be evacuated by the vacuum accumulator. 

2. Fill bellows with Molten sodium - Close all valves except valve D. 

The volume of molten sodium flowing into the bellows is indicated 
by the position indicator. When the bellows maximum capacity is 
reached valve D is closed. 
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Remove trapped air inside receiver tank - The receiver tank is 
connected to the system fill line shown in the Figure 12. To 
evacuate the tank valves B and C remain closed and valve A is opened. 

4. Fill receiver tank with molten sodium - Close valves A and D, actuate 
bellows with nitrogen gas pressure and pump molten sodium into tank. 
To obtain the required sodium purity the sodium is pumped through a 
30 micron and 10 micron filter. 

After the receiver tank is filled to the correct level the intake pipe 
is capped-off and the tank is agitated for approximately 15 minutes. During 
this time, the oxide layer on the inside surface of the receiver will be 
removed by the hot sodium and carried in suspension. After approximately 30 
minutes the hot sodium and oxide impurities are removed from the receiver 
tank leaving an oxide-free surface. The sodium removed from the receiver tank 
can be filtered and reused. After completing steps C and D once again, the 
receiver tank is filled and ready to be mated with the other components. 

6.0 CONCLUSIONS 

1. Sodium has been in extensive use as a heat transfer medium throughout 
the United States and Europe for many years. 

2. The high heat transfer coefficient of sodium and lew vapor pressure 
(less than 1 atmosphere) make it possible to design a simple and 
efficient solar receiver. 

3. The corrosion behavior of stainless steel alloys in hot sodium is 
well studied and very predictable. The low corrosion rates of 
stainless steel make it an excellent material for a long life (15 
years) solar receiver. In addition, stainless steel material is 
relatively inexpensive and readily available. 

4. The handling, filling, filtering, and operating procedures associated 
with hot sodium are very well established. The impurity level in 
the FACC baseline receiver can be held to as low as 10 ppm. 

5. Technical Report number 015 shows that hot sodium is a very safe 
medium when handled properly. 
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